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ABSTRACT 



A semi-empirical method for modelling the loss of electron fluxes in the earth’s mag- 
netosphere was developed. An equation for the integral-energy omnidirectional electron flux 
as a function of time and magnetic field strength was derived from pitch-angle diffusion 
theory. 

This flux equation was the basis for a computer datarfitting program written at the 
Air Force Weapons Laboratory (AFWL) to fit the AFWL Trapped Electron Data Base. The 
program utilised a least-squares fit and incorporated random variations of the characteristic 
exponential loss times about their initial values. An improved table of initial less times was 
compiled for use with the program. 

The derived flux model showed substantial agreement with the empirical data base. 
Representative plots of computed flux over raw data are shown for L-values of 2.3 and 2.4. 



I. INTRODUCTION 



Objective 

Modern technology, both civilian and military, depends greatly on the use of satellite 
systems. It is thus important to the Department of Defense to understand and predict the 
environment in which satellites may operate, in order to determine survivability and vul- 
nerability requirements. One part of that environment is the electron lluxes which may exist 
in the trapping regions of the earth’s magnetosphere. Such (luxes may arise from natural 
sources or from injection by high-altitude nuclear bursts. Because such fluxes have great 
capability for inflicting damage to sensitive satellite components, even at great distances 
from a burst, understanding their behavior is vitally important. 

There are a large number of theoretical treatments in the literature which describe 
various electron loss mechanisms from the earth’s magnetosphere. However, the empirical 
calculations of actual electron losses are often conflicting, and frequently are restricted in 
their areas of coverage. 

Since the cessation of high altitude nuclear testing in 1962, loss studies have been limited 
to electrons injected into the magnetosphere by natural magnetic disturbances such as solar 
storms. Thus, much of the theory which has been developed since 1962 has had restricted 
opportunity for empirical testing. 

The objective of this study is to produce an improved algorithm to calculate losses 
of electrons from the magnetosphere. The algorithm is semi-empirical, because it utilizes 
the Air Force Weapons Laboratory (AFWL) Trapped Electron Data Base (Ref. 37) measure- 
ment of electron fluxes (following the old high altitude nuclear bursts) to determine the 
characteristic loss times of those electrons. 

The loss times thus determined and the algorithm for flux “decay” will then be incor- 
porated by AFWL into the SPECTER computer code (Ref. 7-9) for improved flux calcula- 
tions of satellite environments. 



Scope 

This study is limited to investigation of one loss mechanism, pitch- angle difTusion, using 
a theoretical formalism developed principally by Shulz (Ref. 47,48). The study has also been 
necessarily limited to those regions of space and those energies covered by the satellites from 
which the data base was complied. The loss equations thus developed in this study are only 
implicitly functions of energy, E, and L-coordinate (explained in Chapter II). 

The sheer volume of data and the limited time for this study also limited the testing 
of the algorithm developed to a representative subset of the entire data base. 
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Approach 



A specific solution of the pitch- angle diffusion equation is derived using the method 
proposed by Shulz (Ref. 47). This solution is then fitted to the electron flux data using 
a least-squares fitting program writteu by AFWL, with the coordination of this author. 
The solutions are performed for integral fluxes (above threshold energy of the satellite 
detector) and for constant L-value. The solutions are thus explicitly functions of magnetic 
field strength, 13, but implicitly functions of fi and L . 

Use of the AFWL/NTCTS computer facilities was essential to this study. Well ovpt 72 
hours of actual computer-processing time were used, and over 1000 computer-produced fiux 
plots were examined. 



Presentation 

A brief Background section is presented, which covers some physical relations and 
terminology used in the study of space physics. This section is used to establish uniform 
symbol notation, and to provide a common starting point for the reader unfamiliar with the 
specialized language of space physics. It may be skipped by the reader more familiar with 
the topic. 

Chapter III outlines the theory used to develop the flux equations which were fitted to 
the data. 

Following the theory chapter is a chapter outlining the method of data analysis. 
Descriptions are given of the data base and of the AFWL computer codes used in this study. 

Representative results of the study are presented in Chapter V, in both tabular and 
graphic form. Again, the sheer volume of data makes inclusion of all plots and fitted values 
prohibitive. 

Conclusions and recommendations for improvements for further study are presented in 
Chapter VI. 

The Bibliography includes 65 sources and provides a comprehensive summary of existing 
literature of relevant topics about radiation trapped in the magnetosphere and of pitch-angle 
diffusion. 



D. BACKGROUND 



The equation of motion of a charged particle in magnetic, electric, and external gravita- 
tional fields is: 

f-r 



m~ = mtj-rq 



(f x») + ,S 



where 



r = f(f, ?o, t/ 0 ); — = v 
m = particle mass, 

5 =s gravitational constant, 
q = charge, 

B — magnetic field, 

E = total electric field (Ref. 11:23-24.) 



In a uniform magnetic field, a charged particle will move in helical fashion along and 
around the field lines with a cyclotron radius, p C} (also called gyroradius or Larmor radius) 
about its “guiding center”. The particle’s motion may thus be separated into its rotation 
around the field line and the motion of its center of rotation, or guiding center, along the 
field line (Ref. 11:24-25). 

The period of cyclotron rotation, r e , is defined by (Ref. 39:5): 

2 xm 2 xp r 

Tc ~ qB 



where m is the relativistic mass: 



TTi = rri() "7 = nu) 



1 



and where is component of velocity perpendicular to B. The cyclotron frequency is 
simply 2x over the period (Ref. 39:5): 



_ 2t qB 

T e m 



The angle between a particle's local velocity, v , and the magnetic field, /i, is its pitch- 
angle, a, which is defined as: 



f v \\ 

a = arccos — 

V ” 

In a uniform static magnetic field, is constant and v±_ is constant, and hence a. is constant. 
The particle will then move with a uniform circular motion around the field line and uniform 
rectilinear motion along the field line, resulting in a helical motion about the line, since 

dp fdt = qv X B = ma, 

and 

maji = q(v X B ) j = 




0 



3 n 

w n 

m3j_ 

a ± 



= 0 , 

= eonttant, 

= 7(” X «)_[_, 

= (qv±_B)/m = constant, which is a constant centripetal acceleration 

(Ref.39:4-7; Ref. 11:23-30). 



The concept of the first adiabatic invariant arises naturally from the gi iding center 
approximation. The first adiabatic invariant, M, (also called the relativisi ic magnetic 
moment) is deGncd by: 

ft 

M = — = constant 

2moB 

where pj_ is the particle's perpendicular momentum in the guiding center approximation. 
The assumptions implicit in calling M invariant are that the spatial variation )f B is small 
compared to p e ajid the time variation of B is small compared to r c (Ref. 39:1 1-23). 



If a particle has a constant velocity along a field line in the guiding center frame of 
reference, then 



sin 3 or(/i) sin 2 a,* 

— ----- = — — = constant 
B(s) B { 



where s is the arc length along the field line and : is any point on the field line. 'This assumes 
that the particle’s kinetic energy remains constant as it follows the field line. If a particle 
enters the Geld at point B t - with pitch angle a t -, its velocity along the field line is 



v\\(s) = v cosa(a) = v • 




sin‘a,-. 



If the Geld is increasing in the direction of ujj, then the “mirror point” of the particle, where 
its parallel velocity is aero, is 

n Bi _ 

sin'a,* sin"a(«) 

and the particle has a local pitch angle of 90° at that mirror point (Ref. 39:34 ~42). 



When the magnetic field has a geometry like that of a dipole field, increasing at the 
poles and decreasing in the midpoint, then the particle is in the so-called “magi etic bottle”, 
and is trapped between the mirror points. There is some minimum B value between the 
mirror points which is called Bo, the equatorial field strength. Using the gi iding center 
approximation and approximating the earth’s magnetic field by a dipole fie! 1, a particle 
trapped in the earth’s field has three distinct motions. It circles rapidly about ;. field line, it 
bounces along the line between mirror points, and it drifts slowly in longitude around the 
earth. /Ml three motions take place with different speeds, so they are distinct (Ref. 11:25; 
Ref. 29:34-65). 



From the mirror point definition, and from the velocity equations, we see that 



Note that the particle's mirror point depends only on its “injection” point and its 
injection pitch angle, not on its initial energy or velocity. A mirror point is a consequence 
of the field alone, and all particles injected at the same point with the same pitch angle will 
mirror at the same point on the field line. Of course, this is only true if no external forces 
are acting. 

The parallel velocity, up, of the particle will be the maximun at B 0} where B(a) is the 
minimum. Thus, a trapped particle spends most of its time near the mirror po nts, and the 
least amount of its time transiting the equatorial Geld regions (Ref. 29:34-44). 

The bounce period, r^, oF a trapped particle is generally much greater than its cyclotron 
period, r CJ and is defined by: 

> da 2 da 

= 2 / — 7“7 = - / 

*'•"» v !lW y/\-(B(s)/B m ) 

where a mi and a m2 are the mirror points on the field line. 



As a particle in a dipole field like the earth's bounces along the field line, it also drifts 
perpendicularly to the Geld line, due to external forces, field gradients, and fie’d curvature, 
as well as other effects such as time-dependent field changes. This drift is slew compared 
to the bounce period. As the particle bounces and drifts, it traces out a surface* between its 
mirror points and around the earth, called a drift shell (Ref. 39:9-19). 



The concepts of the second and third adiabatic invariants arise from the above behavior. 
If the forces acting on the particle remain almost constant over its bounce period, the second 
adiabatic invariant , J , is defined by 



/ = 



/ 



P|| * = 




PH da 



where pp is the momentum component parallel to B and da measures arc lcnj;th from the 
equator (Ref. 45:12; Ref. 6:3-31). If the forces acting on a particle remain almost constant 
over its drift period, the third adiabatic invariant, which is the magnetic flu* enclosed by 
a drift shell, is defined by 

4 > = j) Aq • da 

where Ao is the magnetic vector potential and where the integration is performed over a 
curve, a, which lies in the drift shell. The third invariant, 4>, is defined and computed for 
the drift shell of the guiding center, with constant field, and not for a drift shell which a 
particle may physically trace out under short-term conditions (Ref. 45:12; Ref 39:76-79). 



The guiding drift shell of a particle may also be referred to as an “invariant surface,” 
which is composed of field lines which end at the mirror points. The three adiaba. ic invariants 
uniquely define an invariant surface (Ref. 5:3-31 to 3-32). 

The more common set of parameters used to define a particle’s position in the mag- 
netosphere (or to define an invariant surface) is the B-L coordinate system. The B parameter 
is magnetic field strength, and L is the Mcllwain L parameter defined by L = tq j Re where 
Re is the earth’s radius (^ 6371 km) (Ref. 6:3-33) and (Ref. 39:53). 

The value of L , therefore, is equal to the distance, in earth radii, of U c equatorial 
point on a Geld line in a dipole field. If the field is not symmetric (not a pe: feet dipole), 
the invariant surface is not so well defined. However, for most Geld lines in .he trapping 



region of the earth’s magnetosphere, L varies by less than one percent along the line, so 
the B-L system is adequate. The less-than-exact symmetry of the earth’s dipole field and 
outside forces may cause particle drift to deviate from a perfect azimuthal course. This 
means that its L value will vary in the course of a drift period. This variation of L values 
is not significant below L == 3, but may be so above that value. The average lvalue of all 
intersecting invariant surfaces is called an L-skell. This L - value defines a set of surfaces along 
field lines, which may end at diliering mirror points. The //-shells are considered to intersect 
the earth’s surface, even though mirror points do not extend through the atmo; phere to the 
surface (Ref. 0:3-33 to 3-35; Ref. 39:53). 

The various mechanisms which operate to cause trapped particles to change //-shell or 
to be lost from the trapping region will necessarily violate one or more of the adiabatic 
invariants. For example, pitch-angle diffusion violates M or 7, or both; and radial diffusion 
violates 4> (Ref. 45:48). These concepts are discussed further in Chapter 111. 
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m. THEORY 



While the motion of particles trapped in a magnetic field is well understood, the 
mechanisms of losses from magnetic fields such as the earth’s magnetosphere are less well 
understood. It is generally agreed, however, that pitch-angle diffusion into the loss cone is 
one of the predominant mechanisms for removal of charged particles for mid-range Is- values. 
The loss cone angle is the lower limit of pitch angle for trapped particles. Any particles with 
smaller pitch angles will mirror in the sensible atmosphere and will be lost by atmospheric 
scattering. 

Roberts (Ref. 41) has noted that a pitch-angle scattering mechanism must be extant 
for pitch-angle diffusion to occur. Such a mechanism would necessarily violate one or more 
of the adiabatic invariants. 



It is not the purpose of this study to ascertain the true source mechanhm of pitch- 
angle diffusion. Such a study is beyond the scope of this paper. Indeed, the search for 
this mechanism has been going on for decades. Roberts (Ref. 41) has postulated that the 
“perturbation forces’’ causing such diffusion may result from turbulent ambient electric or 
magnetic* fields or from collisions with other trapped or non-trapped particles. Lyons (Refs. 
20-28) has formulated extensive theory based on resonant interactions of so-calh d “whistler” 
VLF waves with the trapped-particles’ gyrofrequencies: “cyclotron-resonance”. Regardless 
of the actual physical mechanism behind pitch-angle diffusion, its general treatment is 
mathematically the same, and the physical results are the same. (Ref. 41:308) This study 
follows the methods of Roberts (Ref. 41), Shull and Lanzerotti (Ref. 45) and Shulz (Ref. 
47) in assuming that pitch-angle diffusion is an operative process. One of the purposes of 
this study is the formulation of the equations necessary to validate that assumption against 
experimental data, and to perform that experimental validation. 

It should be emphasized that, while pitch-angle diffusion is assumed in this paper to 
be the dominant loss mechanism for L-values and altitudes considered, it is not the only 
loss mechanism. Several investigators have formulated radial diffusion (cros:;-L) theories 
(for example, Walt (Ref. 58) or Tomassian (Ref. 54)), coulombic or collisioi-scattering- 
diffusion theories (for example, Wentworth (Ref. 60)) and multiple diffusion theories, such as 
combinations of pitch-angle, energy, and/or radial diffusion (for example, Walt |Ref. 58) and 
Lyons (Ref. 26)), all of which show some agreement with experimental data. la particular, 
atmospheric scattering is obviously a dominant force at very low L values (Ref. 30). 

The introduction of multiple loss-mechanisms makes explicit solutions of any diffusion 
equation extremely difficult. To simplify the problem and render it amenable to the method 
of Shulz (Ref. 47), this study assumes that “as a rule, radial diffusion enables the radiation 
belts to become populated from an external source (or rearranges particles injected by an 
internal source), while pitch angle diffusion causes particle loss to an atmospheric sink” (Ref. 
45:48). 

Pitch-angle diffusion, while simple in concept, is complicated in detail. Roberts’ (Ref. 
41:307-337) treatment of the general mechanism is particularly descriptive: 



At the magnetic equator, a particle’s pitch angle, ao, is determined by 



x = cos a o 



ft 



+ p±, 2 ) 




(i) 



where Pfj 0 and are the parallel and perpendicluar components of the particle's 
momentum with respect to the Geld line at the equatorial value of Dq. Roberts next defines 



(2) 



Wr. 



W 



-Lo 



_ V 

jT7lo 

_ P±l 

2 mo 



where mo is the rest mass. These are defined since the first adiabatic invariant, Af, is 
proportional to 0 : 

VK, 

3) JW = ^ 

Bo 

and, if Wj| w Wj_ # , the second adiabatic invariant, /, is approximately prcportional to 

Ww 

'bWu 

H) J ~ -T- 



where 



rj? = the particle’s bounce period, 
relativistic mass 



7 = 



rest mass 



m 

T7lo 



A particle’s path may be defined by plotting its values of W± 9 and Wj| 0 as in Figure 1. 




Fig. 1. Loss of a trapped particle by a random walk into the pitch- 

angle loss cone. (Ref. 41:306) 

Note that lines of constant equatorial pitch angle pass through the origin, because 

tan 2 a 0 = — 1 . 

W U 
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Any force which violates either M or /, the first or second adiabatic ir variant, will 
produce “diffusion” of the particle’s pitch angle as shown in Figure 1 as W^ 0 and W(| 0 
change with time. When a particle’s equatorial Ditch angle diffuses to the vah e of the loss 
cone angle, 

( 5 ) x c = cos a c = J\ - -jp- 

Y oioo 

where B \ oo is the value at 100 km altitude, then it is lost by atmospheric scattering. “The 
loss cone serves as a ‘sink’ for particles” undergoing pitch angle diffusion (Ref. 41:307). 

Roberts makes two very important points about pitch-angle diffusion as a oss-mechan- 

ism: 

(1) with no source, the entire radiation-belt would be depleted of particlcr, 

(2) the loss cone approximation is just that, since the atmosphere is not shr .rply defined 
at 100 km. However, if atmospheric pitch-angle scattering is not the primary focus, 
the loss come approximation may be useful. “Naturally, when the loss- one concept 
is used, detailed agreement between theory and experiment cannot be expected in. 
the region near the edge of the loss cone” (Ref. 41:308). 

Both Roberts (Ref. 41) and Shulz and Lanzerotti (Ref. 45) develop the pitch angle 
diffusion equation from the Fokker-Planck equation using a particle distribution function, 

7- 



The Fokker-Planck equation is a formalism which arose from the study of Brownian 
motion, and which is used frequently in transport theory (Ref. 42:308). The characteristic 
Fokker-Planck equation for trapped particles (which ignores radial diffusion) it. 



d l 

dt 



1 0 



1 



IP dE [ 7P ( dt )„ + x T(y) dx 



x T{y) D z 



d l 

dx 



Jl JL 

B + IP dE 



7 P Dee 



dE 



where the first term, subscripted v, represents non-stochastic (mean) energy loss to the 
atmosphere, the second term represents pitch-angle diffusion, and the third ter n represents 
range-straggling (energy diffusion). The term T(y) is defined below, and the Dee and D zx 
terms are the characteristic diffusion coefficients; 7 is relativistic mass ratio (R< f. 45:55-58). 



Roberts uses a distribution of particles in a “tube” of force about a fi< Id line, and 
Shulz and Lanzerotti use a phase- space density distribution function which s essentially 
equivalent to Roberts’. Since this study follows the methods of Shulz (Ref. 4") and Shulz 
and Lanzerotti, a discussion of phase space is necessary (Ref. 45:15-22). 

Any moving particle can be described by specifying its three position coordinates and 
its three canonical momentum components. This completes a six-dimension “nhaso-space” 
in which a particle can travel in time. 

If there exists a system of a large number of particles in phase-space, th< system can 
be described by a six-dimensional distribution function /(/*/» <?»»*) where 

Pi (t = 1,2,3) are canonical momentum components, 

(fi (t = 1, 2,3) are position coordinates, and t is time. 
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Thus, / d^Pd z l f is the number of particles instantaneously occupying 6-D volume 
d*Pd 3 q. According to Liouville’s theorem, the phase-space volume containing the system 
of particles moves incompressibly through phase-space. 

Since P is an awkward quantity to deal with physically, Shulz and Lanzerctti note that 
P transforms to the more familiar p as 

(8) ? = P + (!)A 

where A is the electromagnetic vector potential. 



Hence, since q = ? and the P to p transformation has a unit 

Jacobian. The position-momentum distribution function /(p,r, £) defines the particles oc- 
cupying the 6-D volume d z pd z r at any given t. The pitch angle diffusion equation is given in 
terms of a phase-space distribution function which is numerically equal to the more easily- 
definable position-momentum distribution function (Ref. 45:15-22). At constant energy and 
L-shcll, and under the action of some source S, the diffusion equation in phase- >pace can be 
written 



(') 



2 

dt 



l d_ 
X T(y) dx 



xT{y) D xx 



d 2 

dx 



+ 5 



where T(t/) is the quarter-bounce integral path length function of y = >/T - x 2 where 
x = cos oq . The quarter bounce integral gives the length of the trajectory of a particle (in 
units of Rq t the distance to the equatorial crossing) from the equator to the mirror point. 
The exact definition of T is 

r - /'" 

Jo cos a 



where s is the distance along the field line and a is the local pitch angle (Ref. 10:4029—4030). 



D xx is the bounce-averaged diffusion coefficient as a function of x (Ref. 41 and Ref. 
45). The boundary conditions are that / = 0 at x = x e (some cutoff value) and f{x = 0) is 
finite. 



Shulz and Lanzerotti (Ref. 45:162) rewrite the above as 

($) 



d 2 

dt 



1 d__ 
x dx 



x D xx 



d 2 

dx 



- X -D 



n?) 
” nu) 



Ox 



+ 5 



and further approximate that the second term is negligible for x C 1 since t 2 + y 2 = 1. 
If T(y) « T(l), the second term disappears and the equation is a diffusion equation in 
cylindrical coordinates. (These approximations are not used in this work. This development 
is used to aid in understanding the exact solutions which follow from Shulz Ref. 47)). If 
D xx and 5 are then independent of x, and if f(x y t) == X(x) T(t) } the eigenfur ctions of (8) 
are Bessel functions of order zero. The general solution to (8) would then be 



(») 



/(*. o = /«,(*) + a " w J ° 

n 




where is the steady-state solution 



( 10 ) 
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The /e n are the reroes of Jo (n = 1,2,...) and the a n (t) vary as e~W Tn \ where r n = 
x 2 /(D 2 * 2 ). Thus, the pitch angle distribution (and hence the directional or omnidirectional 
flux) is shown to be the sum of a steady state and higher order eigenmodes (Ref. 45:160-168). 

The steady state can be thought of as the “normal” or “quiet-time” value of flux which 
exists in equilibrium with the source 5. An “injection” of particles would thea result in a 
perturbed distribution function with several eigenmodes, each with a characteristic decay- 
time r n . The higher order modes decay faster, and eventually only the fundamental mode 
would remain, which would decay exponentially to reach the steady state (after “infinite” 
time) (Ref. 41 and 45). There is ample evidence that such a process does indeed occur. 
Roberts (Ref. 41) cites Explorer XV data for the 28 October 1062 Russian explosion as an 
example. Rosen and Sanders (Ref. 43) also note that decay is faster immediately after solar 
magnetic storm activity than during quieter periods. The similarities in Figures 2 and 3 
illustrate the way in which a decay of higher eigenmodes (Figure 2) can approximate the 
temporal evaluation of an actual electron distribution (Figure 3). 

The primary difficulties with the above discussion arise from the approximations that: 

(1) T(y) = T(l) = constant 

(2) x 2 <K. 1 (or y « 1) 

(3) D xx varies little (or not at all) with z. 

Roberts (Ref. 41:31 1) uses the argument that the full spiral path (and hence T(y)) varies 
only by a factor of 1.4 as z varies from 0 to 0.9 and by 1.9 as z varies from 0 to 1. Shulz 
and Lanzerotti (Ref. 45:163) use the simplifying assumptions that T(y) = T(l) and D xs is 
independent of x. Both use the assumption that x 2 <£ 1. 

The present study requires that x be allowed to vary from 0 to z c (up to near 1) in 
order to adequately examine the AFWL trapped electron data base (Ref. 37). In addition, 
if x is not small, the approximation that T(y) = T(i) is a poor one, since y is not close to 
1 . 



The function T(y) can be shown, within 0.577S to be approximated as (Rrf. 10:4030) 



(u) 



T(y) « r(0)-|r(0)-T(l)]!/ 3 ^ 



where T(0) = 1 J - 
•15:19). 



(ln(2 + V3))/(2 




1.3802 and T(l) = (*V5)/6 0.7405 (Ref. 



Obviously, if x approaches 1, y approaches 0 and the r(j/) = T{ 1) approximation 
becomes invalid. 



The three limiting approximations are removed in the treatment of Shulz (Ref. 47) by 
introduction of a new “canonical” variable z such that: 

(12) * - Z( V ) = f z 1 TV) dx’ - / ./ 7V) dy' 

J 0 j y 

and a corresponding diffusion coefficient 



(13) 



D:z — [r TV)] 2 D JX 
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Figure 2. Postulated Decay cf Higher Flux Eigenmodes to the Steady Stele by Pitch- 
Angle DiiTusion (Ref. 41:313). 




Figure 3. Temporal evolution of 1.0 Mev Electrons from the Russian 3 Burst fRef. 
41:310). 
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Using (11) in (12) and performing the integration: 



2 




T( 0 )-(T( 0 )-r(l)]/ /4 






(14) 



2 « 



« / y [ r (0)y' - (T(0) - T(l)]/ /4 drj' 

« J r (0)/-^(T(0)-T(l)]/ ,/4 |i 

r(o)(i-» 5 ) 't[T(o)-r(i)][i-y ll / < ] 
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where Z( 0) = 16/35 and Z( 1) = 0 (Ref. 46:5213). 

Using (7), (12) and (13), the diffusion equation becomes: 



(15) 



*2 - L\ D °J. 

Bt dz[ ” d: 



+ s . 



Shall (Ref. 47:6) states that (15) is a “canonical” diffusion equation in that there is no 
Jacobian factor which fails to commute with D xx as was the case in (7). 

The assumptions which now must be made are: 

(1) D xz is some “suitably simple function of z n (Ref. 47:0) so that exact eigenfunctions 
0 rt (c) may be specified, 

(2) 3, the distributed natural environment source, is independent of z . 

Assumption (1) is the limiting assumption, since if the function of z is very :omplex, the 
diffusion equation becomes extremely difficult to solve. Assumption (2) is both simple and 
reasonable, however. The source of the natural environment must be close tc constant in 
order to be the “driver” of quiet-time equilibrium. It must also be distributed airly closely 
to the steady-state quiet time distribution for the same reason. 



Shull (Ref. 47:7) states that even if D xx is not of a functional form tc yield exact 
eigenfunctions, there may be a D xz which resembles D xx “closely enough” and for which 
exact eigenfunctions j/ n (z) (resembling £„(*)) are known. 

The following derivation of the exact eigenfunctions and the corresponding omnidirec- 
tional differential-energy flux arises directly from Shull (Ref. 47). 



The first assumption by Shull is that 




where a is some number, not necessarily an integer, less than 2, and where D xz is the value 
of D zx at some z c < 10/35 where / vanishes. Roberts (1069) also uses this asuimption. 

This form of D xx allows a basis set of orthogonal eigenfunctions to be si own for the 
interval of interest: 0 < z < z e . Following the notation above, we seek seme D xx and 



K 



eigenfunctions 0 n (xr) which “closely resemble” the true functions and which sa isfy 



(17) 

for 

with the substitution 
where 



d 

dz 



D zs 



djM 

dz 



+ ^n5„(~) = 0 



C„ - 

?*(-) = z a w{p: ( ) 

1 — (7 

a — — - — 



and 



/? = 



\ ~ a 

A n ~e 



D„,e 

then (17) becomes Bessel’s equation. The exact eigenfunctions of (17) are given by 
(18) 






whore k utx are the nth zeroes of the Bessel function J u of the order v, where 

cr — 1 

" = 2~—lr' 



The Tj n given by (18) arc normalized so that 

(18) a [ 9n 9m dz = 

J 0 



The eigenvalues of (18) are given by 

I") 



Shulz and Boucher (Ref. 48:6) state that, since particles are not lost at z ~ 0 (but 
rather at z = z e ) } the diffusion current must go to zero at z = 0. A diffusion current at 
2=0 would imply diffusion of pitch angles (and hence mirror points) into the equator, 
which is not reasonable. The condition of no diffusion current at z = 0 co -responds to 
lirn z _ 0 D ZJ g' n (z) = 0. 

Shulz and Boucher (Ref. 48) use the series expansion of J u from Abramowic r , and Stegun 
(Ref. 1 ) to show the limit as z approaches zero. This author shows the limit directly for the 
special case eigenfunction {J~x/o) which is used in this report. The derivatior is shown in 
Appendix A. 
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In order to arrive at the differential-energy omnidirectional flux, Shulz (Ref. 47) assumes 
that f(z ) t)-* foo{z) as t-+ oo; that is, the solution distribution function has some steady-state 
value. This also assumes that S is constant in time , but not in x or z. The distribution 
is necessary in order to expand 7 as a scries of eigenfunctions in the flux equation. To arrive 
at an expression for /oo we start by setting Of/Ot to zero in (15) and integrating twice, 
which gives: 

( 20 ) ?«(*) = r 7T- f 5V') dz" a*. 

J x L'x'z* Jn 

Expanding 5 as a series of orthogonal eigenfunctions gives: 

(20a) $(*") = 

a —0 

and using orthogonality (18a) with (20a) gives: 

( 20b ) a m = [’S{J") rj n (z”') dz 1 ". 

Substituting (20a) into (20) gives 

( 20c ) fco( z )= f 7 T—f a " ') dz" dz'. 

Jx Ds,.Jo ^ 

Rearranging the order of summation and integration in (20c) gives: 

( 20d) ?«>(*) = £ a n f TT— f d„(z")dz" dz 1 . 

n-0 Jx Dx ' x ‘ J ° 

Now integrating (17a) twice gives: 

5n(*) = J t Q~f 0 X " 5n(^') dz” dJ 



or since is a constant, 

( 20e ) r i r dz » d j. 

X n Jx U x tx>JQ 
Substitution of (20e) into (20d) gives: 



( 20f ) 



? M- Vo-IsM 

I oo l"') — ^ \ 



n — 0 



Substitution of (20b) into (20f) gives: 

(21) u[z) = f; hQ f e s(z’) s i n {j) d^ 

n-0 X " 70 

where S^/) is positive (hence a source). 

Shulz and Lanzerotti (Ref. 45:39-40) point out that the directional flux j{E\r) = 
P~ /(; , ||>Pj_i'0 w here / is the same distribution function as in this paper. Sh ilz (Ref. 47) 
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points ont that j most be integrated over all angles in the unit sphere in p-sprce to obtain 
the omnidirectional flux J: 



m 

Since 

and 

it follows that 

and that 
Now since 



/•cos a„ 

J = 4xp 2 / / d(cosa). 

Jo 



sin a sm“ao 



B 



y~ = sin' a 



Bo 






cos 



d(cos a) = 



a= f~H£) 

-iffl/floMy 2 ) 



\ /t^wbo) 



cos 



" = .-|(.-' s )=‘ -5« 



if cos o = cos a*, y 2 = y c 2 ; if cos a = 0, y 2 = #o/£ and (22) becomes 



=wf 



-}(D/Do) Mf) 

liB'/B) yj~\ - v *(D/Uo) 
which, upon combining terms and reversing the limits (-): 



(23) 



= 2 Tp 2 / 



(a,/g) (njjhuji f) 

».* y/l - uHDJbT) 



Integrating (23) by parts 



/ 



it dv = ut; — 




where n — J 



(D/ Bp) <j (y 2 ) 




and since 

df _ £/ gz 

tf(jr) c?z <?(y 2 ) 



to 



■will give (since the definite integral evaluates to zero, because f is 0 at y c by definition): 

(W) 



J = " 2 "’ = /„, -f - K t ) 5; 5P) ^ 



But, from (12), dz = — y T(y) dy. So 



_ -y r(y) d, J _ -y r (y) rf y _ -Hy) 

<%-’) <*(*/) ‘-’y d 'j 2 



Hence, (24) ultimately becomes: 
(25) 



= - Up2 L v^Kl) t(,,) f* 'W. 



Shulz then expands 7 ^ a weighted series of eigenfunctions: 

oo 

(2«) f = f 00 (z)+J2A n (E,L-,t)y n (z). 



n — 0 



This allows the flux / in (25) to be written: 
(27) 



J = 2x p 2 



MB'(B) I TpT 



oo j r r> \ 

- X) ME, L;t)J' t y 1 - y 2 (|j r(y) ?'„(*) <y 2 



The following two expressions are derived in Appendix B. 

(28) <» f /-\ r*» 

/'<»(*) = _ 

i-O A " 



(20) 



f'o c(*)= r~SlnV)dV), 

«-0 X " y ° 

«-«FK«®r '•-.(-(j)" ) 






The actual form of the flux equation awaits the choice of c r, which also determines 
the dependence of D 2 j on z. Shulz uses a = 0 for purposes of comparison of several 
computational methods. This choice is also attractive because it makes D xx independent 
of z (see 17(b)) which was a major assumption to begin with: 

(30) D zs = D„ f . 

This choice is made in the present study for the above reasons and also for the sake of 
simplicity in computations. 
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With (7 = 0, the Bessel functios within ei geo functions become order —1/1, which can 
be simply written as: 



(31) 



J-i MR) 




cos (R), 



(32) 



JL ir JR)^-yJ^s\n(R) 



and the zeroes n n occur at (2 n 4* l)(?r/2). 



cos (i?) 

*R-\/2/(xR) 



This siraplifios (32), since the second term 'will always be zero if «/ , _i/2(*'*) ’ s evaluated 
at k.„ = (2 n + IKjt/2). 

Equation (32) becomes 






?r(2n + l)(?r/2) 



sin 



((2n + I);;) 



(33) 



J '- 1 («n) W ^2n + 1 )(t/2) ^ ^ 



n-f-1 



The integral with respect to z 1 within (28) would thus bo evaluated as: 

r *’ \fiTzc \! z ' I z c \j 



i; 



•/*<) 



COS 



~ \j ^ ( lT + i )( 7 / 2 ) »in((2n+l)5) 

= r VV Z cos((2n+ 1 )(t/2z c ) ^ 

Jo -sin((2n + 1)( jt/2)) 

_ \/2/z c sin((2n+ 1 )(t/2)(z/ Z ,)) 

((2n + 1 )(.t/2; c ))(— l) n+l 



0 2n+ i) £-/) dz ' 



(34) 



s/zjz, sin((2n+ 1 )(t/2)) 



((2n+l)(T/2^))(-ir +1 
The X n in (28) would be evaluated as (from (19)): 



(35) 



X n = (2n+l)’^yD„.. 



The flux equation (27), using (31), (33), (34), and (35), will then simplify to the following 
for (7 = 0: 



(30) 



J 



(T riBtjm 

n. 



■5y/|2 /.t,) «in ((2n-H)(ir/2|) 
(2oH-l)l l r/2^H-ir*‘ 



(2n+ l) 2 (ir/2z e f D, x 






s/W' 



\J ir(2n+U(ir/2) 



(- 1 )" 





”■(2**+ \)(x/2){z/: r ) 





o 

r(2n+ l)(x/2)(z/z e ) 





cos ((2n + l)(;r/2)(z/z e )) 



x((2n+ 1 )(t/2)(z/z c )) 2 v /-_ )T 2_ 



)U/*c) 






r oo Wfl 

- 

U-o •'»«’ 



O./B) 



vw^y 



\/2/(x(2n + f)(x/2)) (-1)" 



_1_/ z 

* V ^ 




x(2n+ l)(x/2)(z/z c ) 



COS 



(<—>;£) 






2 z c y Y x(2t» + l)(x/2 )(z/z c ) 

cos((2n+ l)(x/2)(z/z e ) ) 

x((2n + l)(x/2)(z/z e ))' \/2/(x(2n + l)(z/ 



(2n+l)(x/2) 



o>) 



Equation (36) reduces to the following upon simplification and cancellation of terms: 
(37) 

J 



■ Hi) rw ^ ) 




Ml-!/ 
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The first term in Equation (37) represents a steady state value of the flux, and the second 
term represents a sum of eigenmodes. Each A n (Ref, 45:102) has the form a n (/’)e“*^ r "^^. 
As n increases, the r n * s become shorter; i.e., the higher eignmodes decay fader, so that 
eventually only the fundamental mode is left, which decays exponentially. At t — oo, the 
fundamental mode has reached zero, and only the flux remaining is the steady state flux, 
which corresponds to the natural environment background. 

The remaining substitutions in (37) are z Cl z, T(t/), and y e : 



V/'V*-( W 

( 33 ) - = Z( V) - 1(1 - y 2 )T( 0) - £(7(0) - T{l)](l - y* 11 ' 4 ’), 

Zc = Z (y e ) - 1(1 - y e 2 )T(0) - £(7(0) - 7(L)l(l - y e '"'") 

(as shown in Equation (14)) T(y) is given by Equation (11). 

The variable of integration in (37) is taken as y rather than y 2 for rimplicity of 
integration. 

The preceeding derivation has explicitly followed Shulz (1081). However, Shulz did 
not show how the differential-energy flux related to the integral-energy flu> . Shulz and 
Lanzerotti (Ref. 45:163) state that the integral flux “will then scale as" the differential flux, 
without derivation. Since the AFWL trapped electron data base (Ref. 37) exists primarily as 
omnidirectional fluxes integrated over broad energy bandpass, the relation must be derived. 

It will be shown in the following and Lanzcrotti (Ref. 45) are indeed correct: The 
differential-energy flux equation (37) will scale as the integral flux, and the energy dependence 
will be imbedded in the constants for the source and the eigenmodes. 

From Roederer (Ref. 39:86), the integral energy omnidirectional flux can be written: 

(39) J>e= [ J d.E 

Je 

where J is the differential-energy flux, as in the previous equations, and />£ it* the integral 
energy flux above threshold energy E. 

A principal assumption will be that energy and pitch angle are both independent 
variables of the flux. Thus, (26) would really look like 

< 39 /(*,£)» JcoSm /«,,(-) + f ]A n {E,L-,t)r Jn (:) 

n — »0 



and hence in (25) d~j jdz would look like 



(39 b) 



£? 

0 



Z = ?«,(£) ~L,(=) + £ A n {E,L-t) ?'„(-)• 



n*=* 0 



Using (27), (39a) and (39b), Equation (39) becomes, at a given //-shell: 
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( 40 ) 

J>E = ~ 2 *f E P 2 l<^(E)dE f''l'^z)dx-2xf^ p- dE An L {E) tf n (z) dz 

where /L©,! 1 ) aa( ^ (/^(z) represent nil of the angular-dependent functions of x, r, or z shown 
in (25) through (27). 

Using the previously stated assumption (Ref. 45:152-103) that A ni (E) are of the form 
a n (E)e ~ t ^ Tn ^'^ , 40) becomes: 

J>E = -2xf p 2 f COl (E)dE f f 007 {z)dz 

* E J x i 

(41) P 2 ar,(Ey-‘'' T - iE)) dE [ ' Tj' n (z)dx 

„-0 JE J »> 

at a given L-shell. 



To “remove” the energy dependence of the decay-time terms, the expectation value 
of is computed, assuming that the second term in (41) represents a iistribution- 

function of the energies. 



(42) 






J”p ; q w (g)g-/^ g » dE 
f~ P *a n (E)dE 



Then (41) can be written, using (42): 

/ >£ = 2*[- J b ptf^WdE^j'' f’ 00i (z)dx 
(43) -2rf ;[| g p'a^dE^r'^J^ f n (z)dz. 

Now since it is assumed that the energy and angular shape functions are separable, 
the f ool (E) must include only those parts of (28) which have energy dependence. The only 
terms in (28) which could have any energy dependence are S and X„. 

Shulz (Ref. 47:23) notes that, while - f ^(z) resembles go{z) in functional form, /^(z) 
will coincide exactly with go{z) only if 5 is directly proportional to go(z). Fcr simplicity, 
the assumption is made in this study that 5 is constant over the interval 0 < z < z c . In 
physical terms, this really corresponds to a steady-state isotropic source. It is known that 
the natural environment is not isotropic over cither space or time, but the variations are 
generally small compared to the variations from an injection by a nuclear evei t. 

Now if S and X n are the only terms in (28) which may have energy depe idence, then 
? 00l (&) in (*43) would correspond to 

3(g) 

„_0 *"(£)’ 

and since the only portion of the X n which may have any energy dependence from (10)) is 
D ZZe) then should be proportional to 

5(E) 

bzxSE)' 
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where D KMe {E) may be a constant in E. There may, of course, be some other functional 
energy-dependence included in f ocl (E) which is yet-undetermined. 

Equation (43) is generally of the same functional form as equation (37) if the terms in 
square brackets in (43) correspond to the combined constant terms in (37) as demonstrated 
above and as shown in Appendix D. Equation (43) also bears out Shulz and Lanz^rotti's (Ref. 
45:163) prediction that the integral flux scales with the deferential flux; i.e., the “pitch-angle 
shape” dependence is invariant for different in.1- or integral-energy fluxes. 

Since it is impractical to evaluate an infinite number of terms in (37), only the terms 
up to n = 3 will be shown. Shulz and Lanzerotti (Ref. 45:163) state that the I igher modes 
vanish for n > 2. This in fact means that the higher modes decay so rapidly that modes 
higher that n = 3 should not be seen except at very early times. The Anal form of Equation 
(43), taken to n = 3, is shown in Appendix C. This is the form of the flux equation used in 
the data analysis, except that fewer terms were used than are shown in Apperdix C. 

The general solution for the differential-energy flux (which, from (43), is directly propor- 
tional to the integral-energy flux), using Bessel functions of arbitrary order and with cr- 
dependence explicitly shown, is given in Appendix D. 



IV. ANALYSIS OF DATA 



The AFWL Trapped Electron Data Base consists of a set of computer tapes of satellite 
counting data which hare been assembled from historical archives (such as the National 
Space Sciences Data Center, NSSDC) (Ref. 37). The data cover injections of electrons from 
nine high-altitude nuclear detonations shown in Table l (Ref. 6:5-2). 

The Data Base has been organized by Pfitzer into a coherent set of tapes which 
present, by satellite, net electron omnidirectional integral-energy (above threshold) fluxes as 
a function of time, B, and L-shell. The computed errors for each data point, in units of flux, 
are also given. The background which was subtracted to give net flux is also given (Ref. 
37:11, 188-189). 

Additional data have been collected by AFWL from plots published in early literature. 
These plots have been photographically enlarged, and the data points digitized and placed 
in computer files with formats similar to the Pfitzer data (Ref. 12:670-671, Ref. 31:046-468, 
Ref. 38:637-638). 

The satellites which provided the data, and their orbital parameters, are shown in Tabic 
II (Ref. 44:41-45). 

The analysis of the data base required development of three computer programs by 
personnel at AFWL, in coordination with this author: 

(1) “Program DTABASE" to read raw data, organize it for processing, and output it 
into uniformly formatted files (Ref. 15). 

(2) A program to plot the raw data and to plot flux curves generated by the fitting 
program. 

(3) “Program Elcctrofit”, a fitting program to take data points from DTVBASE at a 
given L-shell and energy and to fit them with the theoretical model leveloped in 
this study. It also generates flux curves using the fitted functions at predetermined 
times to compare with the raw data (Ref. 34). 

The first two programs were written with limited input from the author. The third was 
written under close and extensive coordination with this author, as discussed in this chapter. 
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The DTABASE program, written by Miss Cherise Jarrett of Computer Sciences 
Corporation, reads the raw data files, in various formats created by Pfitzer (Ref. 37), satel- 
lite by satellite, and stores them in a three-dimensional array. The points are then ordered 
(within each energy group) by increasing time since burst. To generate a file of points at a 
specified L-vniue, each data point /rvalue is compared to the desired L and is recorded in 
the output file if a match occurs. Additional points are “created'* by linear interpolation of 
two successive data points if they fall on both sides of the desired L and if the data points 
differ in time by no more than a specified value. The interpolation is performed in L } B } 
time, flux, and error. The output file can be restricted to a speciGed time “window” and to 
a specified range of B values, if desired; however, the data were analyzed over the full range 
of B and over all times which existed in the Data Base. The program has tie additional 
capability to * //-average” data points at the same L if their individual B valu< s are within 
a user specified limit of each other and they are within the required time difference. This 
procedure has the effect of eliminating “double” points, and was used for the majority of 
this study with a limit of one percent. Another capability of the program is to interpolate 
to a specified //-value in a manner analogous to the /^-interpolation described above; this 
capability was unused in this study (Itef. 15:1-2). 
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Table j. !Iigh-allitude nuclear detonations. 
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TABLE D 



Satellites and Their Orbital Parameters 



Name 


Launch 

Date 


Failure 

Date 


Decay 

Date 


Period 

(rain.) 


Perigee 

(mi.) 


Apogee 

(mi.) 


Inch 


Alouette 1 


9-28-62 


unk. 


— 


105.4 


020 


033 


80.5 


Alpha 
Upsilon 1 


9-1-62 


— 


10-26-64 


94.4 


189 


4H 


82.8 


Explorer 15 


10-27-62 


2-9-63 


— 


312.0 


194 


107 ")0 


18.0 


Injun 1 


6-29-61 


3-6-63 


— 


103.8 


534 


631 


67.0 


Traac 


11-15-61 


7-62 


— 


105.0 


562 


720 


32.4 


Telstar 


7-10-62 


2-21-63 


— 


157.8 


593 


35(3 


44.8 


Explorer 4 


7-26-58 


10-6-58 


10-23-59 


110.1 


163 


1372 


50.1 



(Ref. 44:41-15) 
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Id order to determine the optimum time difference to use for interpolation across the 
desired L-value, the DTABASE program was modified under direction of this author. The 
program now presents a summary table of the raw data points for each satellite, showing the 
frequency of occurance of time differences between successive data points from 0.1 minute 
to 5.0 minutes, by tenths of a minutes. The cumulative percentage of points with a specified 
time difference or less is also shown. These summaries are presented for the bursts and 
satellites examined as Tables III through XIII (Ref. 15:67-70). Perusal of these tables showed 
that the optimum time difference for interpolation was 2.2 minutes, primarily to include the 
maximum number of points in the Starfish Telstar file. This time difference wa.*; used for all 
data examined in this study. 



A final capability of DTABASE is to add data points to a specified Rvalue file by 
rounding points in L to the desired value. This procedure assumes that poim-s which are 
closer to the desired L than some specified limit may be considered to occur at the desired 
L. Within the errors in satellite measurements, this procedure is reasonable, and in fact 
appears (by data comparison) to have been used by Roberts (Ref. 41:310), although he does 
not so explicitly state. This capability exists in the program, but was not used n this study 
(Ref. 15). 



The plotting program was developed by Mr. John Burgio of AFWL to generate all the 
plots of raw data and fitted curves shown in this report. The program plots raw data points 
beginning at a specified start time since burst (in days) and covering the time period in a 
specified time ‘window* (shown in Figure headings as “TW"). The L-value, equatorial B- 
value, and B-cutoff values are also shown. If the B-range of points is restricted, the minimum 
and maximim values of B are shown. The equatorial B value is computed from L by the 
relation 

Bo = 



0.308 



1 ? * 

The B-cutoff value was initially computed by using the dipole formula (Ref. 30:55) 



tfe = , / 1 = \/(fio/£cu«). 

V LV4-(-VO 

However, it was later determined that this approximation might be less than optimum; so 
a table of L-values versus B cut was used and was linearly interpolated across the tabulated 
L-values to find the B cut value. This table (Table XIV) was prepared by AFWL, from the 
48-term Jensen-Cain model of the magnetosphere (Ref. 16). It was found that use of Table 
XIV produced B cut (or y eut or z cut ) values more indicative of the global cutoff values. This 
table was used thereafter in both the plotting program and the fitting program (discussed 
below). This procedure is reasonable in view of the use of the same Jensen-Cain model to 
produce the B-L coordinates for the Pfitzer Data Base (Ref. 37). 



The plotting program plots flux as a function of z = \/l^-(Bo7B). ft also has the 
capability to plot flux data points versus time over all values of L and B, or flux versus 
energy at specified B or L. Next to the satellite name on each plot which has fitted curves 
is a number in square bracket? which corresponds to the serial number of the run of the 
fitting program which prodeuced the curves. 



This author’s input to the plotting program was limited to supplying header information 
and specifying layout and format. 

The third program used in data analysis is “Program Electrofit," written by Mr. Ilarry- 
Murphy of AFWL with the close coordination of this author (Ref. 34). Electro, fit is a least- 
squares fitting program which uses the functions developed by this autlmr (Appendix C) 



to fit data points which are output from DTABASE in order to compute the characteristic 
decay times (r„) and linear amplitude constants (a n ). These decay times and constants are 
determined at the specific energy and lvalue which are input with the data. 



28 



Table ITT 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE < S ) : TELSTAR 

FOR BURST: STARFISH 



Tire 


BIN 


NO. DELTA Tires 


CUMULATIVE 


(MIN) 


IN BIN 


PERCENTAGE 


0 0 


0 1 


607 


1 4 


O. 1 


0 2 


3 


1 4 


0 2 


0 3 


7 


1 5 


0 3 


0. 4 


33 


1. 5 


0. 4 


0. 5 


0 


1. 5 


0. 5 


0 6 


22 


1. 6 


0 6 


0 7 


9 


1. 6 


0 7 


0 B 


58 


1. 7 


0. B 


0. 9 


2 


1 7 


0 9 


1. 0 


219 


2. 3 


1. 0 


1 . 1 


29 


2 3 


1 . 1 


1 2 


4 


2. 3 


1 . 2 


1. 3 


56 


2. 5 


1. 3 


1. 4 


12 


2. 5 


1 . 4 


1 . 5 


15 


2 5 


1. 5 


1 . 6 


0 


2 5 


1 . 6 


1. 7 


35 


2 6 


1 . 7 


1. B 


11 


2 6 


1 . B 


1 9 


19 


2. 7 


1 . 9 


2. 0 


33421 


Bl. 4 


2 0 


2 1 


1 78 


61 B 


2 1 


2. 2 


26 


Bl 9 


2 2 


2 3 


B 


81. 9 


2. 3 


2. 4 


5 


Bl 9 


2 4 


2 5 


1 


81. 9 


2. 5 


2. 6 


5 


Bl 9 


2. 6 


2. 7 


0 


81 9 


2. 7 


2. 8 


2 


81 9 


2. B 


2. 9 


2 


Bl 9 


2. 9 


3. 0 


50 


82 0 


3 0 


3. 1 


6 


82. 1 


3 1 


3 2 


0 


02 1 


3. 2 


3. 3 


3 


82. 1 


3. 3 


3. 4 


1 


82 1 


3 4 


3. 5 


1 


B2 1 


3. 5 


3. 6 


0 


82. 1 


3. 6 


3 7 


2 


82. 1 


3. 7 


3. B 


1 


82. 1 


3. B 


3. 9 


5 


82 1 


3. 9 


4. 0 


1451 


B5 5 


4. 0 


4. 1 


5 


B5. 5 


4 1 


4 . 2 


3 


85 5 


4 2 


4 3 


0 


B5 5 


4. 3 


4. 4 


0 


B5 5 


4. 4 


4. 5 


0 


85 5 


4. 5 


4. 6 


2 


B5 5 


4 6 


4. 7 


6 


85 5 


4. 7 


4 B 


‘2 


B5 5 


4. B 


4. 9 


2 


B5 5 


4. 9 


5. 0 


23 


65 6 


5. 


0<~ 


61 15 


100 0 
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Tnble IV 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE ( S > : EXPLORER* 

FOR BURST: ARCUS 1.2 .3 



Tire 


BIN 


NO DELTA TIMES 


CUMULATIVE 


(MIN) 


IN BIN 


PERCENT AGE 


0 0 


0 1 


442 


63 


7 


0 1 


0 2 


133 


82 


9 


0 2 


0. 3 


46 


89 


5 


0. 3 


0 4 


10 


90. 


9 


0 4 


0 5 


8 


92 


1 


0. 5 


0 6 


0 


92 


1 


0 6 


0 7 


1 


92 


2 


0. 7 


0. 8 


1 


92 


4 


0 8 


0. 9 


0 


92 


4 


0. 9 


1 . 0 


0 


92 


4 


1 0 


1. 1 


0 


92 


4 


1 . 1 


1. 2 


0 


92 


4 


1 2 


1. 3 


0 


92 


4 


1. 3 


1 . 4 


0 


92 


4 


1 . 4 


1 5 


0 


92 


4 


1 . 5 


1 6 


0 


92 


4 


1 6 


1. 7 


0 


92 


4 


1. 7 


1 8 


0 


92. 


4 


1. 8 


1 9 


0 


92. 


4 


1. 9 


2 0 


0 


92. 


4 


2 0 


2. 1 


0 


92. 


4 


2 1 


2. 2 


0 


92. 


4 


2 2 


2. 3 


0 


92 


4 


2 3 


2. 4 


1 


92 


5 


2. 4 


2. 5 


0 


92. 


5 


2 5 


2 6 


0 


92 


5 


2 6 


2 7 


0 


92. 


5 


2 7 


2. 8 


1 


92 


7 


2 8 


2 9 


*1 


92. 


8 


2 9 


3. 0 


0 


92 


8 


3. 0 


3. 1 


1 


92 


9 


3. 1 


3 2 


1 


93 


1 


3. 2 


3. 3 


0 


93 


1 


3. 3 


3 4 


0 


93 


1 


3 4 


3 5 


0 


93 


1 


3 5 


3. 6 


0 


93 


1 


3 6 


3 7 


0 


93. 


1 


3 7 


3 8 


0 


93. 


1 


3 8 


3 9 


0 


93 


1 


3. 9 


4. 0 


0 


93 


1 


4 0 


4 1 


0 


93 


1 


4 1 


4 2 


0 


93. 


1 


4 2 


4 3 


0 


93 


1 


4 3 


.4 4 


0 


93. 


1 


4 4 


4 5 


0 


93 


1 


4 5 


4 6 


0 


93. 


1 


4 6 


4 7 


0 


93 


1 


4 7 


4 8 


0 


93 


1 


4 8 


4 9 


0 


93 


1 


4. 9 


5. 0 


0 


93 


1 


5. 


H 

V 

o 


48 


100 


0 



Table V 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE < S ) : ALLOUETTE1 

FOR BURST: RUSSIAN 1 



Tire 


BIN 


NO. DELTA Tires 


CUMULATIVE 




< MI N > 


IN BIN 


PERCENT ACE 


0 


0 


0. 1 


0 


0 


0 


0, 


1 


0 2 


3413 


72 


7 


0 


2 


0. 3 


0 


72 


7 


0. 


3 


0 4 


572 


84 


0 


0. 


4 


0. 5 


180 


88 


7 


0. 


5 


0. 6 


15 


89 


0 


0 


6 


0 7 


105 


91 


2 


0. 


7 


0 0 


0 


91 


2 


0, 


6 


0 9 


46 


92. 


2 


0. 


9 


1 . 0 


46 


93 


2 


1 


0 


1. 1 


1 


93 


2 


1 . 


1 


1 . 2 


5 


93 


3 


1. 


2 


1. 3 


0 


93. 


3 


1. 


3 


1. 4 


15 


93. 


6 


1 . 


4 


1. 5 


4 


93 


7 


1 . 


5 


1 . 6 


3 


93. 


8 


1 . 


6 


1 . 7 


5 


93 


9 


1 . 


7 


1. 8 


0 


93. 


9 


1 . 


6 


1. 9 


12 


94 


1 


1 


9 


2 0 


1 


94. 


2 


2. 


0 


2. 1 


1 


94 


2 


2 


1 


2. 2 


2 


94 


2 


2 


2 


2 3 


0 


94 


2 


2 


3 


2. 4 


2 


94 


3 


2 


4 


2 5 


1 


94. 


3 


2. 


5 


2. 6 


0 


94 


3 


2. 


6 


2. 7 


3 


94 


4 


2. 


7 


2. 8 


0 


94 


4 


2 


0 


2. 9 


0 


94 


4 


2 


9 


3. 0 


0 


94. 


4 


3 


0 


3. 1 


1 


94. 


4 


3 


1 


3. 2 


0 


94 


4 


3. 


2 


3. 3 


0 


94. 


4 


3 


3 


3. 4 


3 


94. 


4 


3. 


4 


3. 5 


0 


94. 


4 


3 


5 


3. 6 


0 


94. 


4 


3. 


6 


3. 7 


0 


94 


4 


3. 


7 


3. 0 


0 


94 


4 


3. 


0 


3 9 


1 


94 


5 


3 


9 


4 0 


1 


94 


5 


4. 


0 


4. 1 


0 


94 


5 


4. 


1 


4. 2 


3 


94 


5 


4 


2 


4. 3 


0 


94 


5 


4 


. 3 


4. 4 


1 


94 


6 


4 


4 


4 5 


0 


94 


6 


4 


. 5 


4. 6 


0 


94 


6 


4 


. 6 


4. 7 


0 


94 


6 


4 


. 7 


4. 0 


0 


94 


6 


4 


0 


4. 9 


0 


94 


6 


4 


. 9 


5 0 


0 


94 


6 




5. 


0<C*= 


255 


100 


0 



Table VI 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE ( S ) : EXPLORER 15 

FOR BURST RUSSIAN 1 



TIME BIN 
(MIN ) 



NO. DELTA TIMES 
IN BIN 



CUMULATIVE 

PERCENTAGE 



0 0 


0 1 


0 


0. 0 


0. 1 


0 2 


0 


0 0 


0 2 


0 3 


0 


0 0 


0. 3 


0. 4 


2 


1 7 


0. 4 


0 5 


12 


12 0 


0. 5 


0 6 


0 


12 0 


0 6 


0. 7 


3 


16 2 


0. 7 


0 8 


6 


21 4 


0. 8 


0 9 


1 


22 2 


0 9 


1 0 


0 


22 2 


1 . 0 


1 1 


13 


33. 3 


1 . 1 


1 2 


8 


40 2 


1 2 


1 3 


0 


40 2 


1 . 3 


1 4 


1 


41 0 


1 4 


1 5 


0 


41 0 


1 . 5 


1 6 


0 


41 0 


1 . 6 


1 7 


11 


50 4 


1 7 


1 8 


3 


53 0 


1. 8 


1. 9 


0 


53 0 


1 9 


2 0 


0 


53 0 


2 0 


2 1 


0 


53 0 


2. 1 


2 2 


3 


55 6 


2 2 


2. 3 


4 


59 0 


2 3 


2. 4 


0 


59 0 


2 4 


2 5 


0 


59 0 


2 5 


2. 6 


0 


39. 0 


2. 6 


2 7 


0 


59. 0 


2 7 


2 8 


0 


59 0 


2 8 


2 9 


‘l 


59. 8 


2 9 


3 0 


0 


59 8 


3 0 


3 1 


2 


61 5 


3. 1 


3 2 


0 


61. 5 


3 2 


3 3 


0 


61 5 


3 3 


3 4 


3 


64. 1 


3 4 


3 5 


7 


70 1 


3 5 


3 6 


0 


70. 1 


3 6 


3 7 


3 


72 6 


3. 7 


3 8 


0 


72 6 


3 8 


3 9 


0 


72. 6 


3 9 


4 0 


1 


73 5 


4 0 


4 1 


0 


73 5 


4 1 


4 2 


0 


73 5 


4 2 


4 3 


0 


73 5 


4 3 


4 4 


0 


73 5 


4 4 


4 3 


0 


73 5 


4 5 


4 6 


1 


74 4 


4 6 


4 7 


0 


74 4 


4 7 


4 B 


0 


74 4 


4 8 


4 9 


0 


74 4 


4. 9 


5. 0 


0 


74 4 


5. 


0< = 


30 


100 0 



Table VTI 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE < S ) : TELSTAR 

FOR BURST; RUSSIAN 1 



TIME 


BIN 


NO. DELTA TIMES 


CUMULATIVE 


(MIN) 


IN BIN 


PERCENTAGE 


0 0 


0. 1 


0 


0 0 


0. 1 


0 2 


0 


0 0 


0 2 


0 3 


0 


0 0 


0. 3 


0 4 


0 


0. 0 


0. 4 


0. 5 


1 


0 0 


0 5 


0 6 


0 


0 0 


0 6 


0. 7 


0 


0 0 


0 7 


0. 0 


0 


0. 0 


0 0 


0. 9 


0 


0. 0 


0. 9 


1 0 


7 


0. 2 


1 . 0 


1 1 


1 


0. 2 


1 1 


1. 2 


0 


0. 2 


1. 2 


1. 3 


0 


0. 2 


1 . 3 


1 . 4 


0 


0. 2 


1 . 4 


1. 5 


0 


0. 2 


1. 5 


1 . 6 


0 


0. 2 


1 6 


1 7 


0 


0. 2 


1 . 7 


1 9 


1 


0 2 


1. 8 


1 9 


2 


0 3 


1. 9 


2. 0 


3637 


86. 6 


2. 0 


2. 1 


3 


86 6 


2 1 


2 2 


3 


86 7 


2. 2 


2. 3 


0 


86 7 


2. 3 


2 4 


0 


86 7 


2. 4 


2 b 


0 


86 7 


2. 5 


2. 6 


0 


86. 7 


2. 6 


2. 7 


0 


86. 7 


2 7 


2. 0 


0 


86 7 


2. 8 


2. 9 


o 


86 7 


2. 9 


3. 0 


6 


86 0 


3 0 


3. 1 


0 


86. 0 


3 1 


3. 2 


0 


86 0 


3 2 


3. 3 


0 


86. 0 


3. 3 


3. 4 


0 


86 8 


3. 4 


3 5 


0 


86 8 


3 5 


3. 6 


0 


86 0 


3 6 


3 7 


0 


86 0 


3. 7 


3. 0 


0 


86 0 


3. 9 


3 9 


1 


86 9 


3. 9 


4 0 


134 


90 0 


4. 0 


4 1 


0 


90 0 


4. 1 


4. 2 


0 


90 0 


4. 2 


4 3 


0 


90 0 


4 3 


4. 4 


2 


90 1 


4 4 


4. b 


0 


90. 1 


4 5 


4. 6 


0 


90. 1 


4. 6 


4 7 


0 


90 1 


4 7 


4. 8 


0 


90. 1 


4 8 


4 . 9 


0 


90 1 


4. 9 


5. 0 


5 


90 2 


5. 


O 

A 

II 


413 


100. 0 






Table VTJI 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE ( S ) : ALLOUETTE1 

FOR BURST RUSSIAN 2 

TIhE BIN NO. DELTA TirtS CUMULATIVE 

(MIN) IN BIN PERCENTAGE 



0. 


0 


0 


1 


45 


1 


1 


0 


1 


0 


2 


3107 


78 


9 


0 


2 


0 


3 


0 


78 


9 


0 


3 


0 


4 


409 


89 


1 


0 


4 


0 


5 


91 


91 


4 


0 


5 


0 


6 


7 


91 


5 


0 


6 


0 


7 


52 


92 


8 


0 


7 


0 


8 


0 


92 


8 


0 


8 


0 


9 


29 


93 


6 


0 


9 


1 


0 


25 


94 


2 


1 


0 


1 


1 


3 


94 


3 


1 


1 


1 


2 


13 


94 


6 


1 


2 


1 . 


3 


0 


94 


6 


1. 


3 


1 . 


4 


7 


94 


8 


1. 


4 


1 


5 


1 


94 


8 


1 


5 


1 


6 


0 


94 


8 


1 


6 


1 


7 


6 


94 


9 


1. 


7 


1 


8 


0 


94. 


9 


1 


8 


1 


9 


4 


95 


0 


1 


9 


2. 


0 


1 


95 


1 


2. 


0 


2 


1 


0 


95. 


1 


2 


1 


2. 


2 


4 


95. 


2 


2. 


2 


2 


3 


0 


95. 


2 


2 


3 


2 


4 


1 


95 


2 


2 


4 


2 


5 


3 


95 


3 


2 


5 


2 


6 


0 


95 


3 


2 


6 


2. 


7 


0 


95 


3 


2 


7 


2 


8 


0 


95 


3 


2 


8 


2 


9 


0 


95 


3 


2 


9 


3. 


0 


0 


95 


3 


3. 


0 


3 


1 


0 


95 


3 


3 


1 


3 


2 


0 


95 


3 


3. 


2 


3. 


3 


0 


95 


3 


3 


3 


3 


4 


0 


95 


3 


3 


4 


3 


5 


0 


95 


3 


3 


5 


3 


6 


0 


95 


3 


3 


6 


3 


7 


0 


95 


3 


3 


. 7 


3 


8 


0 


95 


3 


3 


8 


3 


9 


0 


95 


3 


3 


9 


4 


0 


0 


95 


3 


4 


0 


4 


. 1 


0 


95 


3 


4 


1 


4 


2 


0 


95 


3 


4 


2 


4 


3 


0 


95 


3 


4 


3 


4 


4 


0 


95 


3 


4 


. 4 


4 


5 


0 


95 


3 


4 


5 


4 


6 


0 


95. 


3 


4 


6 


4 


7 


0 


95 


3 


4 


7 


4 


8 


0 


95. 


3 


4 


8 


4 


9 


0 


95 


3 


4 


9 


5 


0 


0 


95 


3 




5 


O 

A 




189 


100 


0 



?A 



Table DC 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE ( S ) : EXPLORER 15 

FOR BURST; RUSSIAN 2 



TIME BIN NO. DELTA TIMES CUMULATIVE 

(MIN) IN BIN PERCENTAGE 



0, 0 


0 1 


0 


0 0 


0 1 


0 2 


0 


0 0 


0 2 


0 3 


0 


0 0 


0. 3 


0. 4 


71 


2. 9 


0. 4 


0. 5 


264 


13. 8 


0. 5 


0 6 


344 


28 0 


0 6 


0 7 


142 


33 9 


0 7 


0. 8 


77 


37. 0 


0 0 


0 9 


31 


30. 3 


0 9 


1. 0 


83 


41 7 


1 0 


1. 1 


65 


44. 4 


1. 1 


1. 2 


60 


47 2 


1. 2 


1. 3 


64 


49 9 


1 3 


1 . 4 


63 


52. 5 


1 . 4 


1. 5 


115 


57 2 


1 5 


1. 6 


75 


60 3 


1 6 


1 . 7 


125 


65 4 


1 7 


1. 8 


36 


66 9 


1 8 


1 9 


1 1 


67 4 


1 . 9 


2. 0 


33 


68 7 


2. 0 


2. 1 


12 


69 2 


2 1 


2 2 


29 


70 4 


2. 2 


2. 3 


32 


71 8 


2. 3 


2. 4 


17 


72 5 


2. 4 


2. 5 


42 


74 2 


2. 5 


2 6 


14 


74 0 


2. 6 


2 7 


26 


75. 0 


2 „ 7 


2. 8 


11 


76. 3 


2. 8 


2. 9 


16 


76. 9 


2. 9 


3. 0 


10 


77. 4 


3. 0 


3 1 


10 


77 8 


3. 1 


3 2 


10 


78 2 


3 2 


3 3 


7 


78 5 


3 3 


3 4 


8 


78 0 


3. 4 


3 5 


22 


79 7 


3. 5 


3. 6 


23 


80. 7 


3. 6 


3. 7 


45 


82 5 


3. 7 


3. 8 


14 


83 1 


3. 8 


3. 9 


18 


83 8 


3 9 


4. 0 


27 


84. 9 


4 0 


4 1 


15 


85. 6 


4 1 


4. 2 


10 


86. 0 


4. 2 


4. 3 


6 


86 2 


4. 3 


4. 4 


10 


86 6 


4 . 4 


4. 5 


4 


86 0 


4. 5 


4 6 


15 


87 4 


4. 6 


4. 7 


10 


87. 8 


4. 7 


4 8 


4 


88 0 


4 8 


4. 9 


14 


88 6 


4. 9 


5 0 


0 


88. 9 


5. 


0<= 


269 


10O. 0 
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Tabic X 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL ITE < S ) : TELSTAR 

FOR BURST. RUSSIAN 2 

TICC BIN NO. DELTA TIfES CUMULATIVE 

(MIN) IN BIN PERCENTAGE 



0 0 


0 1 


0 


0. 0 


0 1 


0 2 


0 


0 0 


0 2 


0 3 


0 


0. 0 


0 3 


0. 4 


0 


0 0 


0. 4 


0 5 


0 


0 0 


0 5 


0 6 


0 


0 0 


0 6 


0 7 


0 


0 0 


0. 7 


0 B 


0 


0. 0 


0 B 


0. 9 


0 


0. 0 


0 9 


1 0 


2 


0 1 


1. 0 


1 1 


0 


0 1 


1 . 1 


1. 2 


0 


0 1 


1 2 


1 3 


0 


0. 1 


1 . 3 


1 4 


0 


0 1 


1 4 


1 5 


0 


0 1 


1 . 5 


1 . 6 


0 


0 1 


1 . 6 


1 7 


0 


0 1 


1 7 


1 8 


1 


0 1 


1 B 


1. 9 


0 


0. 1 


1 . 9 


2. 0 


2244 


85. 6 


2. 0 


2. 1 


49 


87. 5 


2 1 


2. 2 


1 


B7. 5 


2. 2 


2. 3 


0 


B7. 5 


2 3 


2 4 


0 


67 5 


2. 4 


2. 5 


0 


B7. 5 


2. 5 


2 6 


0 


87. 5 


2. 6 


2 7 


0 


87. 5 


2. 7 


2. B 


0 


87. 5 


2. B 


2 9 


*0 


87. 5 


2. 9 


3 0 


3 


87 6 


3. 0 


3 1 


0 


87 6 


3. 1 


3 2 


0 


87. 6 


3 2 


3 3 


0 


87 6 


3 3 


3. 4 


0 


87 6 


3 4 


3 5 


0 


87. 6 


3. 5 


3 6 


0 


87. 6 


3. 6 


3. 7 


0 


87. 6 


3. 7 


3 . 8 


0 


87 6 


3. B 


3. 9 


2 


87. 7 


3. 9 


4 0 


100 


91 5 


4. 0 


4. 1 


2 


91 6 


4. 1 


4 2 


0 


91 6 


4 2 


4 3 


0 


91 6 


4 3 


4 4 


0 


91 6 


4 4 


4 5 


0 


91. 6 


4 5 


4 6 


0 


91 6 


4. 6 


4 7 


0 


91 6 


4 7 


4. 8 


0 


91 6 


4 B 


4. 9 


0 


91 6 


4 9 


5 0 


0 


91. 6 


5. 


OO 


221 


100 0 



Tnble XI 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE ( S ) : ALLOLCT TE1 

FOR BURST: RUSSIAN 3 



T IhE 


BIN 


NO. DELTA TIMES 


CUMULATIVE 


(MIN) 


IN BIN 


PERCENTAGE 


0 0 


0 1 


123 


0 8 


0 1 


0 2 


10901 


74 0 


0 2 


0 3 


45 


74 3 


0. 3 


0. 4 


1595 


85 0 


0. 4 


0. 5 


491 


80 3 


0. 5 


0. 6 


55 


88 7 


0. 6 


0. 7 


289 


90 6 


0. 7 


0. 8 


11 


90 7 


0. 8 


0 9 


128 


91 5 


0. 9 


1 0 


74 


92. 0 


1 0 


1. 1 


8 


92 1 


1 1 


1. 2 


55 


92 5 


1. 2 


1. 3 


3 


92 5 


1. 3 


1. 4 


62 


92. 9 


1 . 4 


1. 5 


28 


93 1 


1 . 5 


1. 6 


9 


93. 1 


1. 6 


1. 7 


10 


93. 3 


1. 7 


1. 8 


0 


93 3 


1. 8 


1. 9 


10 


93 4 


1. 9 


2. 0 


15 


93. 5 


2. 0 


2. 1 


2 


93 5 


2. 1 


2. 2 


10 


93 6 


2. 2 


2 3 


0 


93 6 


2 3 


2. 4 


13 


93. 7 


2. 4 


2. 5 


7 


93 7 


2. 5 


2. 6 


1 


93. 7 


2. 6 


2 . 7 


7 


93. 8 


2. 7 


2. 8 


0 


93 8 


2. 8 


2. 9 


*4 


93. 8 


2. 9 


3. 0 


2 


93 8 


3 0 


3 1 


2 


93 8 


3 1 


3 2 


5 


93 8 


3. 2 


3. 3 


0 


93 8 


3. 3 


3. 4 


2 


93 9 


3. 4 


3. 5 


0 


93 9 


3. 5 


3. 6 


1 


93 9 


3. 6 


3. 7 


1 


93 9 


3. 7 


3. 8 


0 


93 9 


3 8 


3 9 


1 


93. 9 


3. 9 


4. 0 


0 


93 9 


4 0 


4. 1 


0 


93. 9 


4. 1 


4. 2 


0 


93. 9 


4 . 2 


4. 3 


0 


93. 9 


4. 3 


4. 4 


0 


93. 9 


4. 4 


4. 5 


0 


93. 9 


4. 5 


4. 6 


1 


93 9 


4. 6 


4. 7 


0 


93 9 


4. .7 


4. 8 


* 0 


93 9 


4 8 


4. 9 


0 


93 9 


4. 9 


5. 0 


0 


93. 9 


5. 


O 

A 

N 


911 


100. 0 



Table Xn 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL ITE < S ) : EXPLORER 1 5 

FOR BURST: RUSSIAN 3 



T I 


rE 


BIN 


NO. DELTA TIhtS 


CUMULAT 


IVE 




<MIN> 




IN BIN 


PERCENT 


ACE 


0 


0 


0 


1 


46 


0 


3 


0 


1 


0 


2 


0 


0 


3 


0. 


2 


0 


3 


77 


0 


0 


0. 


3 


0 


4 


500 


4 


0 


0. 


4 


0 


5 


1654 


16 


0 


0. 


5 


0. 


6 


2144 


30 


5 


0 


6 


0 


7 


1629 


41 


5 


0 


7 


0 


0 


970 


40 


0 


0 


0 


0 


9 


274 


49. 


9 


0 


9 


1. 


0 


194 


51 


2 


1 


0 


1 


1 


226 


52 


7 


1 


1 


1. 


2 


210 


54. 


1 


1. 


2 


1. 


3 


293 


56. 


1 


1. 


3 


1 


4 


273 


50 


0 


1. 


4 


1 


5 


474 


61 


2 


1. 


5 


1 


6 


442 


64 


1 


1. 


6 


1 


7 


610 


60 


3 


1 


7 


1. 


0 


313 


70 


4 


1. 


0 


1. 


9 


194 


71 


7 


1. 


9 


2 


0 


131 


72 


6 


2 


0 


2. 


1 


02 


73 


2 


2 


1 


2 


2 


92 


73 


0 


2 


2 


2 


3 


74 


74 


3 


2. 


3 


2 


4 


54 


74 


7 


2 


4 


2. 


5 


45 


75 


0 


2. 


5 


2 


6 


35 


75 


2 


2 


6 


2. 


7 


53 


75 


6 


2. 


7 


2. 


0 


34 


75 


0 


2 


e 


2. 


9 


37 


76 


0 


2 


9 


3 


0 


40 


76. 


3 


3. 


0 


3 


1 


29 


76 


5 


3 


l 


3. 


2 


47 


76 


0 


3 


2 


3 


3 


52 


77 


2 


3 


3 


3 


4 


66 


77 


6 


3 


4 


3 


5 


53 


70 


0 


3 


5 


3 


6 


79 


70 


5 


3 


. 6 


3 


7 


155 


79 


6 


3 


. 7 


3 


0 


91 


80 


2 


3 


. 0 


3. 


9 


127 


01. 


0 


3 


9 


4 


0 


1 1 1 


01 


0 


4 


0 


4 


. 1 


03 


02 


3 


4 


1 


4 


2 


1 17 


03 


1 


4 


2 


4 


3 


87 


83 


7 


4 


3 


4 


4 


1 10 


04. 


5 


4 


4 


4 


. 5 


73 


05 


0 


4 


5 


4 


. 6 


100 


05. 


6 


4 


6 


4 


. 7 


76 


86 


2 


4 


7 


4 


0 


63 


86 


6 


4 


0 


4 


. 9 


00 


87 


1 


4 


9 


5 


. 0 


59 


87 


5 




5. 


o 

A 

H 




1040 


100 


0 



ns 



Table Xni 



COUNTS OF TIME DIFFERENCES BETWEEN NONZERO RAW DATA POINTS 
FOR ALL CHANNELS OF SATELL I TE ( S ) : TELSTAR 

FOR BURST: RUSSIAN 3 



TIME BIN 
(MIN) 



NO. DELTA TIMES 
IN BIN 



CUMULATIVE 

PERCENTAGE 



0. 0 


0. 1 


5 


0 0 


0 1 


0. 2 


0 


0 0 


0. 2 


0. 3 


0 


0. 0 


0. 3 


0. 4 


0 


0. 0 


0. 4 


0. 5 


0 


0. 0 


0. 5 


0. 6 


0 


0. 0 


0. 6 


0. 7 


0 


0. 0 


0. 7 


0. 8 


0 


0 0 


0 0 


0. 9 


1 


0. 0 


0. 9 


1. 0 


64 


0 3 


1. 0 


1. 1 


4 


0. 3 


1. 1 


1 2 


2 


0 3 


1. 2 


1. 3 


0 


0. 3 


1 3 


1. 4 


0 


0 3 


1. 4 


1. 5 


1 


0. 3 


1. 5 


1. 6 


3 


0. 3 


1. 6 


1. 7 


3 


0. 3 


1. 7 


1. 8 


6 


0. 4 


1 . 8 


1. 9 


15 


0. 4 


1 9 


2. 0 


20107 


82 9 


2. 0 


2. 1 


121 


83 4 


2. 1 


2. 2 


7 


83 5 


2 2 


2. 3 


6 


83 5 


2 3 


2. 4 


1 


83 5 


2. 4 


2. 5 


0 


83 5 


2. 5 


2 6 


0 


83 5 


2. 6 


2„ 7 


0 


83 5 


2 7 


2. 8 


0 


83 5 


2. 8 


2 9 


0 


B3 5 


2. 9 


3. 0 


25 


83 6 


3 0 


3. 1 


1 


83 6 


3 1 


3. 2 


3 


83 6 


3 2 


3. 3 


1 


83 6 


3. 3 


3. 4 


0 


83 6 


3. 4 


3. 5 


0 


83 6 


3. 5 


3. 6 


0 


83 6 


3. 6 


3. 7 


3 


83 6 


3. 7 


3. 8 


2 


83 6 


3. 8 


3. 9 


9 


83 7 


3 9 


4 . 0 


1027 


87 9 


4 0 


4. 1 


5 


87 9 


4. 1 


4 . 2 


1 


87. 9 


4. 2 


4 . 3 


0 


87 9 


4 3 


4. 4 


3 


87 9 


4. 4 


4. 5 


1 


87 9 


4. 5 


4 6 


0 


87 9 


4. 6 


4. 7 


0 


87 9 


4. 7 


4. 8 


' 0 


87 9 


4. 0 


4 9 


1 


87 9 


4 9 


5. 0 


12 


88 0 


5. 


0<=» 


29 29 


100. 0 






Table XIV 

Global Magnetic Cutoff Values 



L- Value 


Bcutof f 


1.09*19 
1.1152 
1.1356 
1.1550 
1.1762 
1.2628 
1 1817 
1.8331 
2.3168 
2.9330 
3.6816 
4.5626 
5.5760 
6.7218 
8.0000 


0.23001 

.22908 

.22182 

.22288 

.22353 

.24017 

.26770 

.27573 

.29762 

.33311 

.37001 

.40274 

.43178 

.45502 

.46951 



Note: The Table may be interpolated for other ZrTaJues (Ref. 16). 
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The overall flow for the program is as follows (see Appendix E): 

(1) Use an initial guess for the fundamental decay time and the theoretical ratios 
of the higher-mode de^y times to compute the starting decay tin es. The first 
guess at fundamental decay time is shown in Table XV as a function of L and 
energy. This table is interpolated in L and E for starting value. 

(2) Numerically compute the integrals of the functions shown in Appemlix C at each 

data point using the initial guessed decay times. The integration s performed 
by a 12-point Gauss/Legcndrc Quadrature at each data point for each function. 
Only the fundamental decay mode and two higher modes were computed. 

(3) Compute the linear amplitude constants using a least-squares matrix-solving 
routine. Compute the standard error of the fit. The individual points are 
weighted by the reciprocal of their fractional errors, which have been normal- 
ized to an average value of one. 

(4) Optimize the initially computed y eu t value by using a “Golden Section Min- 
imim” function. This section attempts to let the data determine the “best" 
cutoff value. If the data are not distributed down toward the |< ss cone, the 
Golden Minimum function will provide an erroneous y cut value Hence, the 
function specifies that the computed |/ <rut must be within ±10 percent of the 
initial input value. If, in the iteration process, the function tries to exceed these 
limits, the initial guess of y CXAt is speciGed as the correct value. 

(5) Repeat the steps (2) and (3) with the optimum y eut value. 

(0) To determine the best combination of decay times and amplitude constants, 
the decay times are changed by a “random-walk” method, either increased 
or decreased. The amount by which the decay times may randomly vary is 
determined by using a cumulative normal probability distribution about the 
initial value. 

The standard deviation used was estimated according to this author’s best guess of 
the accuracy of the tables and curves shown by Stassinopoulos (Ref. 53:31-32, 40 — 14) and 
West (Ref. 63:50-54). The Stassinopoulos tables do not always correlate with each other 
or with the Stassinopoulos curves, and the West curves do not always correlate with the 
Stassinopoulos in the L-regions of overlap. The initial fundamental decay times in Table XV 
were computed by this author by interpolation of the tables and curves in Reference 53 and 
by reading and extrapolation of the curves in Reference 63. The author’s estimates of their 
standard deviations are shown in Table XVI. Of particular note is the very low confidence 
in decay times at high L values (above 2.2) and at high energies (above 2.0 MeV). There 
have been few computations of decay times in these regions, and the results vary widely 
(Refs. 53, 03). 

The effect of using the cumulative normal probability distribution and the above stan- 
dard deviations is to restrict the “random-walk” of the decay times away from their initial 
values. The farther away from the initial “mean value” that the random val le generator 
places the new decay time, the greater the “push” that is created back toward he mean for 
the next walk. The distance away, however, is defined in terms of the standaid deviation, 
so that r values in which there is low confidence may vary by larger amounts than may r 
values in which there is greater confidence. The r’s are also restricted by the co istraint that 
r 0 > t\ > ro. 
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With each change of decay times, steps (2) and (3) are repeated up to 05 times, in an 
attempt to get a “better” fit by reducing the standard error of the fit. 



Initial Estimate of Fundamental Decay Time (days) 
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(Compute*! from lief. r ,lt .'U 1(1 II and Kef. «:i :50 5 - 1 ) 



Percent Standard Deviation of Initial 
Fundamental Decay Times (in Table XV) 



L- Range 

Energy Range 
MeV 


1.02- 

1.20 


1.21- 

1.40 


1.41- 

1.80 


- § 
OC • 


> 

2.20 


E < 0.0 


50 


30 


20 


35 


200 


0.0 < E < 2.0 


30 


20 


20 


75 


100 


E > 2.0 


50 


20 


50 


100 


100 



(estimated from Ref. 53:31-32, 40-44 and Ref. 53:50-54) 



(7) If a better fit is found in less than 00 random walks, and if less than 32 “bettor” 
values hare been found, steps (2), (3) and (0) are repeated up to 32 times, using 
the latest “best” decay times. 

(8) Up to five repeats of steps (4) and (5) may be carried out, but only if 32 “better” 
fits are found. Otherwise, only the first y cut optimization is performed. 

(9) The program will always end at the point where 96 random-walk3 have failed 
to reduce the error of the fit from the previous “best-fit”. The decay times, 
y cu t , and amplitude constants are printed for each “better” fit, with the final 
values being those of the “best” fit. 

(10) After the best fit is found, the program uses the fitted constants and decay 
times to compute fluxes, at specified times, from the equator to the best cutolf 
value. These specified times correspond to default times plus 1/2 ime window 
in ^he program which plots raw data. Thus, the computed fluxer are used to 
generate curves which are plotted over raw data points at approximately the 
same times. 

Variations in the program were tried in order to study the effects of the fitting process. 
The principal modifications were in the random decay-time variations. 

The first modification was to fix the ratios of the r’s to their theoretical values of 
1:1/9:1/25 (Ref. 48:10-19). The effect of this modification is to significantly reduce the 
degrees of freedom available to fit the data. 

The second modification was a compromise between the above and complete freedom to 
random walk. This variation allowed r 0 and r { to vary freely (subject to the usual contraint 
that r 0 > n), but r x and r 2 were fixed at their theoretical ratio of 25:9. The principal reason 
for this attempt was the very short time-span covered by some satellite/burst combinations. 
For example, the Russian 2 Burst covered only 4 1/2 days prior to the Russian 3 injection. 
If the data cover a time-span which is small compared to the fundamental decay constant, 
then the fitted value of that decay time may be suspect. However, if the time-s >an is not so 
small compared to the decay times of the higher modes, one should have more confidence 
in the fitted values of those decay time. Stated another way, the fit may not be sensitive to 
the fundamental, but may still be sensitive to the higher modes. 

In this analysis, not all L-values covered by the satellite in Table II wore examined. 
Over 1000 plots of flux versus x were examined at various L-values for the satellite/burst 
combinations of Tables I and II. In many cases, the satellite orbits were of such high 
inclination (toward polar) that data were not available close to equatorial B- values. Because 
of this, and because of the very large volume of data, those satellite/burst/L combinations 
which provided the most consistent coverage from the equator to the loss cone were chosen 
for study. 

Generally, the data fell into three regions of x-space (with some overlap). Region 1 was 
from the equator to the point where flux dropped off sharply into the loss cone (x = 0 to 
x £5* 0.8 or 0.8). Region II was the area where the flux “turned the corner” from a rather 
“horizontal” curve to a rather “vertical” curve into the loss cone (x ^ 0.7 to 0.9). Region III 
was the loss cone region, where the flux levels dropped sharply to the cutoff vtlue (x ^ .9 
for a typical L-value). 



Representative plots of raw data in the three regions are shown in Figures 4 through 8. 
These plots are typical in that data seldom cover all three regions in a short time-span. Gaps 
in the data appear because of satellite orbital coverage and because of lack of on-satellite 
recorders. When a satellite lacked a recorder, data were collected only as the satellite passed 
above the radio horizon of a ground tracking station. Figure 9 is an eight-hour plot of typical 
satellite coverage over all B-L space. 



INTEGRAL ELECTRON FLUX (ELECTRONS/CM^)/SEC 



BURST : RUSSIAN 2 SAT : ALLOUETTE1 
TIME : 0 (dy) 0 (hr) 0.00 (min) O.OOO(days) 

L : 2.40 Beq : .022 BCUT : 0.511 TW (days) : 0.5 

Bmin (Gauss) : 0.001 Bmax (Gauss) : 1.000 

10 9 
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io 3 

X = (1-Beq/B) 1/2 

Figure 4(a). R,rw Allouette Data Collected Only in Region III (at burst tine). 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : RUSSIAN 2 SAT : ALLOUETTE1 



TIME : 3 (dy) 0 (hr) 0.00 (min) 3.000(days) 

L : 2.50 Beq : .020 BCUT : 0.515 TW (days) : 1.0 
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Figure 1(b). Raw Allouette Data Collected Only in R.egion III (3 days post-burst). 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



10 



10 



10 



10 



10 F 



10 



10 s 



BURST : STARF SH SAT : T ELS’AR 
TIME : 1 (dy) 0 (hr) 0.00 (min) I.OOC(days) 

L : 1.80 3eq : .053 BCUT : 0.470 TW (days) : 1.0 

Bmin (Gauss) : 0.001 Bmax (Gcuss) : 1.000 
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X = (l-3eq/B) 1/2 

Figure 5. Ra-w Telstar Data Collected in Region* 1 and HI (early time). 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : STARFISH SAT : TELSTAR 
TIME : 15 (dy) 0 (hr) 0.00 (min) 15.000(days) 
L : 1.70 Beq : .063 BCUT : 0.460 TW (days) : 5.0 
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BURST : STARFISH SAT : TELSTAR 
TIME : 80 (dy) 0 (hr) 0.00 (min) 80.000(dcys) 
L : 1.70 Beq : .063 BCUT : 0.460 TW (days) : 5.0 

Bmin (Gauss) : 0.001 Bmax (Gauss) : 1.000 
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INTECRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 
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BURST : STARFISH SAT : TELSTAR 
TIME : 100 (dy) 0 (hr) 0.00 (min) lOO.OOO(days) 
L : 1.80 Beq : .053 BCUT : 0.470 TW (days) : 5.0 

Bmin (Gauss) : 0.001 Bmax (Gauss) : 1.000 
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Figure 8(a). Raw Telstar Data in Regions I and II (late time). 
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INTECRAL ELECTRON ("LUX (ELECTRONS/CM 2 )/SEC 



BURST : STARFISH SAT : TELSTAR 
TIME : 100 (dy) 0 (hr) 0.00 (min) lOO.OCO(days) 
L : 2.00 Beq : .038 BCUT : 0.487 TW (days) : 5.0 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



AFWL TRAPPED ELECTRON DATA BASE 
BURST : STARFISH SATELLITE : TELSTAR 
DAY : 237 HOUR : 1 MIN : 0 jf : 98 

CHANNEL NUMBER : 4 TIME WINDOW : 8.00 

10 8 
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io 5 



io 4 



io 3 

012345678 

TIME (Hours) 

Figure 9. Typical Raw Data Coverage for 8 Hours (B-L values not shown). 




Figures 10 (a) through (g) shows raw data from Explorer XV over the time period of the 
Russian 2 and 3 bursts, at times corresponding to those shown in Figure 3, for comparison. 
It is apparent that fewer poiots are plotted from the AFWL data base in F' gurus 10 (a) 
through (g) than were plotted by Roberts in Figure 3 for the 1.9 Me- electrons, even though 
the satellite and time periods of coverage were the same. These discrepancies arc not resolved 
in this study, although the previously alluded-to L-rounding procedure may account for the 
differences. 

The temporal progression of flux after one burst from one satellite (Telstar), at a typical 
L-value of 1.9, is shown in Figures 11 (a) through (h). Again, typical gaps in coverage of 
Regions I, II, and III are readily apparent over the 90 days shown. 

As previously noted, only fundamental and two higher eignemodes (/o,/i,/ 2 ) *re used 
to fit the satellite data (Appendix C). The third eigenmode (/a) is not used because Shulz and 
Lanserotti (Ref. 45:163) suggest that the modes greater than two rapidly vanish. Initially, the 
steady-state function (foo) was used along with the fundamental, but this solution was found 
to be “competing" with the fit for the fundamental, in the sense that the linear constants 
tended to be roughly equal in magnitude and tended to alternate in sign. Shulz (Ref. 47:22- 
23) indicates that the /oo and fo functions should resemble each other in shape. Also, Pfitzer 
(Ref.37) has subtracted the background from all flux data, and the background should 
approximate the steady-state, /oo, solution. For these reasons, the steady-state solution was 
not used in the data fit. 



INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 
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BURST : RUSSIAN 2 SAT : EXPLORER15/EXPLORER15 
TIME : 0 (dy) 4 (hr) 58.07 (min) 0.207(days) 
L : 1.90 Beq : .045 8CUT : 0.479 TW (days) : 1.0 

Bmin (Gauss) : 0.001 Bmax (Gauss) : 1.000 
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Figure 10(u). Raw Explorer 15 Data at Times Corresponding to Figure 3. 



INTEGRAL ELECTRON FLUX (ELECTRONS/ CM^)/SEC 



BURST : RUSSIAN 2 SAT : EXPIORER15/EXPLORER15 
TIME : 2 (dy) 19 (hr) 22.09 (min) . 2.807(days) 
L : 1.90 Beq : .045 BCUT : 0.479 TW (days) : 2.0 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 
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BURST : RUSSIAN 2 SAT : EXPLORER15/EXPLORER15 
TIME : 9 (dy) 19 (hr) 22.09 (min) 9.807(days) 
L : 1.90 Beq : .045 BCUT : 0.479 TW (days) : 2.0 
Bmin (Gauss) : 0.001 Bmax (Gauss) : 1.000 
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Figure 10(c). Raw Explorer 15 Data at Times Corresponding to Figure 3. 






INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : RUSSIAN 2 SAT : EXPIORER15/EXPLORER15 
TIME : 20 (dy) 19 (hr) 22.09 (min) 20.807(days) 
L : 1.90 Beq : .045 BCUT : 0.479 TW (days) : 2.0 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : RUSSIAN 2 SAT : EXPLORER 15/EXPLORER15 
TIME : 29 (dy) 19 (hr) 22.09 (min) 29.807(days) 
L : 1.90 Beq : .045 BCUT : 0.479 TW (days) : 4.0 
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Bmin (Gauss) : 0.001 Bmax (Gauss) : 1.000 
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Figure 10(e). Raw Explorer 15 Data at Times Corresponding to Figure 3. 
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INTEGRAL ELECTRON FLUX (ELFCTR0NS/CM 2 )/SEC 
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Figure 10(f). Raw Explorer 15 Data at Times Corresponding to Figure 3. 
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Figure 10(g). Raw Explorer 15 Data at Times Corresponding to Figure 3. 
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Figure 11(a). Temporal Progression for Raw Telstar Data from 1 to 90 Days Post-Burst. 
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Figure 11(b). Temporal Progression for Raw Telstar Data from 1 to 90 Dnyr. Post-Burst. 
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Figure 11(c). Temporal Progression for Raw Tclstar Data from 1 to 90 Days Post-Burst. 
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Figure 11(d). Temporal Progression for Raw Telstar Data from 1 to 90 Days Post-Burst. 



r.O 



INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : STARFISH SAT : TELSTAR 
TIME : 10 (dy) 0 (hr) 0.00 (min) lO.OOO(days) 
L : 1.90 Beq : .045 BCUT : 0.479 TW (days) : 5.0 




10 



10 ° 



■= 10 



10 



10 



10 



10 



Figure 11(e). Temporal Progression for Raw Telstar Data from 1 to 90 Days Post-Burst. 
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Figure 11(g). Temporal Progression for Raw Telstar Data from 1 to 90 Days Post-Burst. 
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V. RESULTS AND DISCUSSION 



Figures 12 and 13 show the typical temporal progression of raw data at two representative 
L values, L = 2.3 and L = 2.1. Included on those figures are the flux curves calculated from 
the amplitude coefficients and decay times fitted to that data by Electrofit. fis previously 
noted, the curves shown on any given figure may vary slightly in time from the individual 
data points, since the curves represent only the one time-instant corresponding to the start 
time plus one-half the time window of the plot. Figures 12 and 13 represent maximum 
freedom of the fit, since each r is free to vary independently of the other r’s. 

Figure 14 shows the same temporal data progression with different fitted curves. The curves 
in this figure were calculated under the restriction that the ratios of the r’s 'vere fixed at 
the theoretical values of 1:1/9:1/25; hence, only T\ was really varying in a random sencc. 

Figure 15 shows similar data progression and fitted curves where the r’s were allowed to 
freely vary, subject to the constraint that Tq > r t > r 2 . Also, this was an initial run with 500 
tries for convergence, vice 96 tries. In addition, the computed fiux was used in the iteration 
process to eliminate non-phys'eal (negative) solutions. The solution was restric ted to being 
positive at four different times and at five different equatorial pitch angles. Alt lough this is 
an initial run with only a small amount of data used, it indicates that the recommendations 
in Chapter VI should be pursued. 

Table XVII shows the fitted linear amplitude coefficients and exponential decay times (in 
days) at constant Energy and L-shell for a representative subset of the data examined in 
this study. The sheer volume of resultant fits and calculated curves made it impossible to 
include all results in this report. The values shown in Table XVII were computed with the 
r’s allowed to vary freely. 

Figures 12 and 13 are considered representative of the results of the fitting model. It is 
apparent from examination of Figures 12 and 13 that the pitch-angle diffusion theory 
and flux computation method developed in this report are consistent with experimental 
measurements. Higher eigenmodes are readily discernible at the lowest energy in the Figures 
or the first five days, and are not visible at later times. These higher eigenmodes are not 
discernible in the plots at the higher energies; Table XVII shows that these modes are 
present, but at much lower levels, in some cases. The calculated curves show a reasonable 
fit to the data even at late times, consistent with the assumption of exponential decay. 

Examination of Figures 14 and 15 shows that some inconsistencies remain to be resolved in 
the use of this fitting method. The caluclatcd flux curves for the low-energy data points also 
appear to be slightly lower than optimum, while the high-energy data appear? to be fitted 
better by the model. This may be partially explained by the existence of many lower data 
points in the vicinity of the loss cone which tend to “drag down” the overall fit. Any change 
in the cutoff value significantly affects the fit by shifting the “vertical” portion of the curve 
toward or away from those points. However, this cannot account for the entire difference 
of the curves from the data. The energy dependence may be treated less than optimally 
in the choice of functional solution made in this study. If the chosen functional solution is 
not the correct one, it may be nevertheless close enough to correct to give a r 2 asonablc fit 
at high energy, but be significantly in error at lower energy. A different functional solution 
may fit the lower-energy groups better and yet still lit the high-energy elect ons as well. 
Recommendation for the use of other functions is made in Chapter VI. 
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Figure 12(a). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2.3 for 
Russian 3 Burst (r„ freely varying). 
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Figure 12(g). Temporal Progression of Raw Data and Fitted Flux Cur-res at L = 2.3 for 
Russian 3 Durst (r„ freely varying). 



7S 



INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : RUSSIAN 3 SAT : EXPLORER15 [045] 
TIME : 25 (dy) 0 (hr) 0.00 (min) 25.000(days) 
L : 2.30 Beq : .025 BCUT : 0.506 TW (days) : 5.0 




10 



10 



10 



10 



10 ' 



10 



10 



10 



Figure 12(h). Temporal Progression of Rot Data and Fitted Flux Curves at L — 2.3 for 
Russian 3 Burst (r n freely varying). 
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Figure 12(i). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2.3 for 
Russian 3 Burst (r„ freely varying). 
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Figure 13(a). Temporal Progression of Raw Data and Fitted Fins Curves at L 2.-1 for 
Russian 3 Durst (r„ freely varying). 
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Figure 13(b). Temporal Progresmon of Raw Data and Fitted Flux Curves a i L — 2.-1 for 
Russian 3 Burst (r* freely varying). 
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Figure 13(c). Temporal Progression of Raw Data and. Fitted Flux Curves at L — 2.4 for 
Russian 3 BuTst (r„ freely varying). 
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Figure 13(g)* Temporal Progression of Raw Data and Fitted Flux Curves at L — 2A for 
Russian 3 Burst (r n freely varying). 
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Figure 13(h). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2.4 for 
Russian 3 Burst (r n freely varying). 
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Figure 14(c). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2.3 for 
Russian 3 Burst (r n fixed 1:1/9:1/25). 
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Figure 14(d). Temporal Progression of Ra^ Data and Fitted Flux Curves at L — 2.3 for 
Russian 3 Burst (r n fixed 1:1/9:1/25). 
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Figure 14(e). Temporal Progression of Raw Data and Fitted Flux Carves at L = 2.3 for 
Russian 3 Burst (r a fixed 1:1/9:1/25). 
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Figure 14(f). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2.3 for 
Russian 3 Burst (r„ O.xcd 1:1/9:1/25). 
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Figure 15(a). Temporal Progression of Raw Data and Fitted Flux Curves at L = 2.1 for 
Russian 3 Durst (r n freely varying with physical contraints). 
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Figure 15(b). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2.4 Tor 
Russian 3 Burst (r a freely varying with physical contraints). 
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Figure 15(c). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2.4 for 
Russian 3 Burst (r n freely varying with physical contrainls). 
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Figure 15(e). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2 A for 
Russian 3 Burst (r„ freely varying with physical eoatraints). 
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Figure 15(f). Temporal Progression of Raw Data and Fitted Flux Curves at L = 2.4 for 
Russian 3 Burst (r n freely varying with physical cootrainls). 
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Figure 15(g). Temporal Progression of Raw Data and Fitted Flux Carres at L = 2.4 for 
Russian 3 Burst (r„ freely varying with physical eonlraiuts). 
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Figure 15(h). Temporal Progression of Raw Data and Fitted Flux Curves at L = 2A for 
U U 35 inn 3 Burst (r n freely varying with physical contmints). 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : RUSSIAN 3 SAT : EXPLORER15 [568] 
TIME : 8 (dy) 0 (hr) 0.00 (min) 8.000(days) 




0 0.5 1 

X = (1-Beq/B) 1/2 
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Figure lS(i). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2.1 for 
Russian ?, Burst (r n freely varying with physical contrnints). 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM^)/SEC 



BURST : RUSSIAN 3 SAT : EXP10RER15 [568] 
TIME : 10 (dy) 0 (hr) 0.00 (min) lO.OOO(days) 
L : 2.40 Beq : .022 Bcut : 0.302 TW (days) : 5.0 




Figure 15(j). Temporal Progression of Raw Data and Fitted Flux Curves at L — 2 A for 
Russian o Burst (r„ freely varying with physical contraints). 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : RUSSIAN 3 SAT : EXPL0RER15 [568] 
TIME : 35 (dy) 0 (hr) 0.00 (min) 35.000(days) 




0 0.5 1 

X = (1-Beq/B) 1/2 
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Figure 15(k). Temporal Progression of Raw Data and Fitted Flux Curv^ at L = 2.4 for 
Russian 3 Burst (r n freely varying with physical contraiots). 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 



BURST : RUSSIAN 3 SAT : EXPLORER15 [568] 
TIME : 50 (dy) 0 (hr) 0.00 (min) 50.000(days) 
L : 2.40 Beq : .022 Bcut : 0.302 TW (days) : 5.0 




Figure 15(1). Temporal Progression of Raw Data and Fitted Flux Curves at L = 2 A for 
Russian 3 Burst (r n freely varying with physical contraints). 
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INTEGRAL ELECTRON FLUX (ELECTRONS/CM 2 )/SEC 






BURST : RUSSIAN 3 SAT : EXPL0RER15 [568] 
TIME : 60 (dy) 0 (hr) 0.00 (min) 60.000(days) 
L : 2.40 Beq : .022 Bcut : 0.302 TW (days) : 5.0 




X = (1-Beq/B) 1/2 
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Figure 15(m), Temporal Progression of Raw Data and Fitted Flux Curves at L ~ 2.4 for 
Russian 3 Burst (r a freely varying with physical coutraints). 
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Pitted Liucar Amplitude Coefficients and 
Exponential Decay Times at constant L-shell and Energy 
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(First line of pair is iuili:il itcrutiou: second line is result.) 
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(|'ir«t line «f p-.tir is nnliiil iteration; second line is 



The short time-spans covered by most of the data sets in the Trapped Electron Data Base 
have not allowed adequate determination of r*). The values of m presented in Table XVII 
are suspect for this reason. However, more confidence is held in the values of t\ and r 2 
computed by this study In most data sets studied, the time period of the data was long 
enough compared to the values of r t and rw to compute reasonable values. More consistency 
between data sets was found for and r 2 than was found for r 0 . Of course, then r 0 may be 
simply computed by multiplying r t by the theoretical ratio corresponding to the solution. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



Conclusions 

The pitch-angle diffusion theory and flux calculation method developed in this study 
are reasonably consistent with experimental data in most cases. The method si ows promise 
as a way of calculating improved decay times over a broad region of the magnetosphere, 
particularly for higher energy electrons. Additionally, a comprehensive summary table of 
the best available literature values for r 0 has been prepared (Table XV). A complete table 
covering such a broad range of energy and L has not heretofore been compiled. 

Time and computer limitations have precluded examination of the entire data base. 
However, the overall efficiency of the model developed in this study has been demonstrated. It 
is expected that improved predictions of fluxes and improved inputs to the AFWd SPECTER 
Codes (Ref. 7-9) will be available as a result of this work. This will enhance the capability 
for calculation of satellite operational environments, and will improve the abili y to predict 
satellite survivability and vulnerability. 



Recommendations 

The following recommendations for further study are proposed to improve the model 
developed in this study, and to resolve the problem of non-physical solutions for some data 
sots. At the direction of this author, several of these recommendations are presently being 
implemented at the Air Force Weapons Laboratory. 

(1) Incorporate an estimation of the error in the magnetic field value, Z?, of 
individual data points. This corresponds to a “horizontal* error bar on the 
plots such as Figure 12. It is known that errors exist in the computation of the 
D value associated with each data point, and this error would prove significant 
in the region of the loss cone. Improved weighting of the data points with the 
combined errors should improve the fit. 

(2) Utilize the /wounding procedure of Chapter IV to add data points to each Is- 
value considered. This should aid in filling gaps in the data, which necessarily 
will improve the fit. 

(3) Discard the highest eigenraode if the fitted value of its decay time is less that 
0.25 day, and fit the data with only the fundamental and one higher mode. 
A value of 0.25 day is not reasonably fitted within the errors in the Trapped 
Electron Data Base. This mode would be essentially lost within th) first day in 
any case. 

(4) Perform an “energy- scaling* of the data to combine data from different satel- 

lites. This involves assuming some form of energy spectrum, such as 
including an estimate of the error in the assumed value of Eq. The new flux 
would be simply the old flux multiplied by a factor £•*<<)/£• ]f the 

error in the old flux were 50 percent (a value consistent with the Data Base) 
and the error in Eq were 25 percent, and if Eq 1.1 (a fission spectrum), then 
the resultant error in the shifted flux would be less than 1 percent different 
from that of the old flux, when shifting from 660 KeV to 500 KeV. Thus, it 
should be possible to combine data sets and fill the gaps in data to achieve a 
better fit. However, the limitation of the assumed spectrum still exists. 
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(5) Utilize only tate-time data where it is available, and Gt only the fundamental 
mode to it, under the assumption that higher modes will have decayed away. 
Then, with the fundamental mode constrained, go back to early time data and 
attempt to fit higher modes, if any. This recommendation is of limited value 
because only two data sets cover a significant period of time: Starfish/Telstar 
and Russian 3/Bxplorer 5. However, it may provide a better value of r 0 than is 
presently available at some L values. 

(f>) Revise the fitting functions to other Bessel function solutions than the J ~\/2 
solution and compare results. This was performed only on one other functional 
solution, Jo, with only one small data set, and the results showed no improve- 
ment over the present study. However, to be certain that the author has not 
simply chosen the wrong functional form, more functional forms must be in- 
vestigated. 

(7) In the fitting program, constrain the directional flux to physical reality at each 
iteration (i.e., compute flux at each iteration and reject non-physical solutions 
at each iteration, rather than after converging on a final solution). This should 
converge to a physically real solution if it converges at all; however, significant 
increases in computer processing time may accrue. 

(8) Increase the number of tries for a solution from 98 to 200 and reduce the number 
of outer loops from 32 to 10 or 15. This may improve the convergence of the 
least squares error, if non-physical solutions are also rejected at each iteration. 

(9) Compare the background subtracted from the flux with the latest available 
natural environment, to ascertain the correctness of the background subtrac- 
tion. If significant differences are found, the background should be re-added to 
the flux and the natural environment subtracted. 

(10) A suggestion proposed by Professor D.G. Shankland, Air Force Institute of 
Technology, involves a significant modification to the fitting method. This 
method involves an interpolating function, or measure of “roughness” of the 
fit, which may be made arbitrarily smooth. The advantage of this method is 
that no limiting assumptions need be made about the data. The chief disad- 
vantage is that the method is a strictly mathematical method, with no physical 
constraints. Although the method would involve extensive reprogramming, it 
promises worthwhile information about the flux curves, and should be tried. 
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APPENDIX A 



Demonstration of Zero Diffusion Current 
at the Magnetic Equator for Half-Integer 
Bessel Function Solution 



To demonstrate zero diffusion current, one must show that 

liin r _o D zx j' n (z) = 0 

for the <7 = 0 case. We Know that from (31): 

cos ( r) 

and from (16): 
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Hence, the derivative of the eigenfunction: 
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Finally, the quantity which must go to zero is: 
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The term in square brackets is a finite constant for any n, and clearly the sine goes to zero 
as z goes to zero. Hence, the limit is zero. 
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APPENDIX B 



Derivation of Eigenfunction Derivatives and y' n 
To find the derivative of the steady state function, / f3el one starts with (21): 



7oo = E 't f S(S)T Jn (;>)d(S) 

„-o X " Jq 



If the X n are assumed (from (19)) independent of z t then the series can simply be 
differentiated term-by-term, since the definite integral is simply a constant for each value 
of n. 
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To find the derivative of the eigenfunction g n , one starts with (18): 






Using the product rule for differentiation: 



APPENDIX C 



The Specific Solution of the Integral- Energy 
Omnidirectional Flux Equation for tbe Half-Integer 
Bessel Function Solution 

The integral-energy omnidirectional flux Equation (43), which is the energy- integrated 
form of Equation (37), is represented as the steady state term, / plus the sum of exponen- 
tially decaying eigenmodes fa through fa: 

j > = V />o7 = foo + fa + 5\ + /2 + /3 



In this appendix, each term is shown, as derived from Equations (37) and (43). Note 
that the eigenmode sum should be infinite, but only the fundamental mode and three higher 
modes are shown. 



The constants in each term are defined as follows (Ref. 39:55): 
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The “fitting constants” are defined as in Equation (43) or in Appendix D: 
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These fitting constants arc really fuctions of energy and //-shell. 

The ratios of ro, T \ , r 2; and r 3 should be (as computed from the eigen vain** ratio of the 
zeroes of the function): 

ro « 9rj r, = (l/0)r n 
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APPENDIX D 



The General Solution of the Integral-Energy 
Omnidirectional Flux Equation for Bessel Functions 
of Arbitrary Order, as a Function of Sigma 

Substitution of (18), (28), and (29) into (27) completes the general expro;sion for 
differential energy omnidirectional flux, J , with Bessel functions of order 



a — 1 




where o is some number less than two: 
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The analogy between the above equation and Equation (43) is seen if: 
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and if: 
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lienee, (43) is simply the energy-integrated form of (37). 
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